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vors. The numerical analysis in this paper demonstrates that the corrections 
at NNLO are sizable, as expected from earlier work. 
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I. INTRODUCTION 

The generally accepted theory of the strong interac- 
tion, Quantum Chromo Dynamics (QCD), provides in 
principle a way to derive various low-energy hadronic ob- 
servables, such as masses and decay constants. However, 
an analytic ab initio derivation of such properties has not 
yet been possible. An alternative approach called Lattice 
QCD, where the functional integrals are evaluated using 
numerical Monte Carlo techniques on a discretized space- 
time lattice, suggests a way around this problem, but is 
nevertheless not free of difficulties of its own. In partic- 
ular, the properties of low-mass particles are especially 
difficult to calculate as they can propagate over large dis- 
tances on the lattice, giving rise to large nonlocal correla- 
tions which limit the precision obtainable with given com- 
puter resources. As a consequence, most simulations have 
been performed with heavier quark masses than those of 
the physical world. The quark masses used in present 
Lattice QCD simulations typically fulfill m u ^ > m s /8. 
The results then have to be extrapolated down to the 
physical masses of ~ m s /25. 

The preferable way to perform the extrapolation to 
physical quark masses is by means of Chiral Perturba- 
tion Theory (xPT) 0, , which provides a theoretically 
correct description of the low-energy properties of QCD 
within the range of quark masses where xPT calcula- 
tions can be considered accurate enough. In particular, 
the appearance of so-called chiral logarithms renders sim- 
ple polynomial extrapolation insufficient. If the Lattice 
QCD simulations are performed with quark masses in 
the chiral regime, i.e. with quark masses for which xPT 
calculations are reliable, then they can be used to e.g. 
determine the effective low-energy constants (LEC:s) of 
%PT. Nevertheless, in practice it is still difficult to reach 
the chiral regime and therefore many large scale lattice 
simulations of so-called quenched QCD have been per- 



formed, in which the effects of closed (sea) quark loops 
are neglected. Consideration of such loop effects requires 
repeated evaluation of fermion determinants, which is 
computationally extremely expensive. It should be noted 
that first Lattice QCD calculations with sea quark masses 
significantly below m s are now becoming available. 

From a computational point of view, it is however 
much cheaper to vary the valence quark masses rn V ai 
rather than those of the sea quarks, m sea . Therefore it 
is not uncommon to perform simulations where the sea 
quark loops are not neglected, but rather suppressed by 
choosing m sca > m va i , or even n sea 7^ ^val > where n va i and 
n sea denote the number of each quark species in the the- 
ory. Such procedures are referred to as partial quenching 
and lead into a space of unphysical theories. 

A. Partially Quenched Theories 

At first, one might expect that partially quenched 
(PQ) theories could yield only qualitative information 
about QCD. However, since unquenched QCD may be re- 
covered from partially quenched QCD (PQQCD) in the 
limit of equal sea and valence quark masses, it follows 
that QCD and PQQCD are continuously connected by 
variation of sea quark masses. In contrast, this is not true 
for fully quenched QCD. Nevertheless, xPT can be ex- 
tended to the case of quenched QCD, which was done by 
Bernard, Golterman and Sharpe 0, y| after earlier work 
by Morel The extension to the partially quenched 
case was done in Ref. 0. A more extensive discussion 
can be found in the work of Sharpe and Shoresh 0, ■ 

The formulation of Partially Quenched Chiral Pertur- 
bation Theory (PQxPT) is such that the dependence on 
the sea quark masses is explicit, and thus the limit of 
equal sea and valence quark masses can also be consid- 
ered for PQxPT. It follows that x?T is recovered as a 
continuous limit of PQ^PT just as QCD is from PQQCD. 
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In particular, the LEC:s of xPT, which are of physi- 
cal significance, can be obtained directly from those of 
PQxPT. In addition, the ability to vary valence and sea 
quark masses separately allows more information to be 
extracted at fewer values of the sea quark masses. This is 
another reason to consider PQQCD in lattice simulations 
as discussed in detail in Ref. [jj. 

This paper is devoted to the extension of PQxPT to 
next-to-next-to- leading-order (NNLO). The discussions 
already given in Refs. 0, 0, are completed and the 
calculations of the masses and decay constants of the 
charged, or off-diagonal, pseudoscalar mesons are com- 
pleted for the case of n soa = 3. It should be noted that 
indications already exist, by the jqq+ql collaboration [l2T ] 
and the MILC collaboration 13], that the 0(p 6 ) contri- 
butions to these quantities are sizable, and that the inclu- 
sion of such effects can have a significant impact on the 
chiral extrapolations down to the physical quark masses. 
More examples can be found in the proceedings of the 
Lattice 2005 conference. 



tation used for loop integrals and combinations of quark 
masses. Sects. 11111 and II VI contain the analytical NNLO 
expressions for the masses and decay constants of the 
charged pseudoscalar mesons. Sect. presents a discus- 
sion of the checks performed on the analytical and nu- 
merical calculations, a numerical analysis of the results, 
and an elaboration on the extraction of the various LEC:s 
from Lattice QCD calculations is given in Sect. IV El Our 
conclusions are shortly discussed in Sect. I VII 

II. TECHNICAL ASPECTS OF PQxPT 

The technical aspects of PQxPT calculations to NLO 
have been thoroughly covered, for the supersymmetric 
formulation of PQxPT, in Refs. 0,11] ■ The new parts spe- 
cific to the NNLO calculations have been covered briefly 
m Refs. UGSIIIl. 

The aim of this section is to collect 
all this information in one place and describe the NNLO 
aspects in somewhat greater detail. 



B. Partially Quenched xPT at NNLO 

The meson masses and decay constants of the pseu- 
doscalar mesons in PQxPT for n sca = 3 to NLO, one 
loop, or 0(p 4 ) in the momentum expansion, were calcu- 
lated in Refs. 0, Q, @ • First results for these quantities 
at NNLO, two loops, or 0(p 6 ), namely the mass of a 
charged pseudoscalar meson for degenerate sea and de- 
generate valence quark masses, can be found in Ref. @. 
The NNLO expression for the decay constant of a charged 
pseudoscalar meson with two nondegenerate sea quarks 
has been calculated in Ref. 01 ■ The full results for pseu- 
doscalar meson masses and decay constants for n sca = 2 
were given in Ref. |ll| . This paper presents the full non- 
degenerate calculations of the charged pseudoscalar me- 
son mass to NNLO for n sca = 3, along with the calcu- 
lation of the decay constant for three nondegenerate sea 
quarks. 

In general, the NNLO expressions are very long, but it 
is possible to shorten them considerably by introducing a 
specialized notation which can accommodate the compli- 
cations inherent in PQxPT. Furthermore, this notation 
satisfies many different algebraic relations, which allows 
for a systematic simplification to be carried out. In the 
next sections, the technical background for the NNLO 
calculations, together with the above-mentioned notation 
is presented, followed by the results for the NNLO masses 
and decay constants of the charged pseudoscalar mesons, 
and a numerical analysis. All analytical formulas are 
given explicitly, but they can also be downloaded from 
the website fl4j . 

This paper is organized in the following manner: 
Sect. UTI introduces PQxPT and discusses the various as- 
pects relevant for calculations at NNLO. In particular, 
the relations between the various sets of LEC:s are high- 
lighted. This section also includes an overview of all no- 



A. Supersymmetric PQxPT 

Let us first recall a few aspects of xPT, and note in 
passing that introductions to the subject can be found 
e.g. in Ref. 15]. For nf massless quark flavors, QCD has 
a chiral symmetry 

G=SU(n f ) L xSU(n f ) R , (1) 

which is spontaneously broken by the vacuum to its di- 
agonal subgroup H — SU(nf)v, where V = L + R. The 
Goldstone bosons produced by this spontaneous break- 
down live on the coset G/H which is also an SU(rif) 
manifold and can be parameterized in terms of a special 
unitary matrix U. This matrix is conventionally written 
in terms of a traceless Hermitian matrix according to 

U = exp (iV2cj)/F} (2) 

where the constant F with dimension of energy can be 
shown to be related to the decay constant of the pion. 
For three quark flavors (u, d, s), the matrix (j) has a flavor 
structure 

(uu ud us \ 
du dd ds I . (3) 
su sd ss J 

The above considerations do not distinguish between 
valence and sea quarks whereas the construction of a PQ 
theory requires a mechanism which gives different masses 
to sea quarks and valence quarks. This effect may be pro- 
duced in a systematical way by adding to xPT explicit 
sea quarks as well as unphysical bosonic ghost quarks. 
The latter cancel exactly all effects of closed valence 
quark loops due to their different statistics, provided that 
their masses are identical to those of the valence quarks. 
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Supersymmetric PQxPT thus contains a set of fermionic 
as well as bosonic quarks. This leads to a modification 
of the chiral symmetry group as given by Eq. (JJJ , which 
in the PQ theory has the graded structure 

G = SU(n val + n sca |n va i) L x SU(n val + n sca ,\n va ,\) R . (4) 

The precise structure of G is somewhat different as dis- 
cussed in 0, This theory contains n va \ valence and 
n^ai ghost quarks, as well as n sea flavors of sea quarks. 
As the PQ theories include bosonic ghost quarks, they 
are not normal relativistic quantum field theories since 
they violate the spin-statistics theorem. However, under 
the assumption that the low-energy structure of such a 
theory can be described similarly to the case of normal 
QCD, one arrives at an effective low-energy theory in 
terms of a matrix U, according to 



U = exp 



(iV2$/F 



(5) 



where the matrix $ now has a more complicated flavor 
structure than </> in Eq. (0 because of the different types 
of quarks present. In terms of a sub-matrix notation for 
the flavor structure 



u a u b u a dp u a Sb 
'In'lb = | d a u b d a dp d a s b 
s a iib s a db s a Sb 



(6) 



the matrix $ becomes 



$ = 



qvqv 



qsqv 



\ qsqv 



qvqs 



qsqs 



qBqs 



qvqs 



qsqs 



qBqs 



(7) 
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where the labels V, S and B stand for valence, sea and 
bosonic quarks, respectively. The size of each sub-matrix 
depends on the exact number of quark flavors used. 

The quarks qy, qs and their respective antiquarks are 
fermions, while the quarks q R and their antiquarks are 
bosons. Thus a combination of a fermionic quark and a 
bosonic quark yield a fermionic (anticommuting) field, 
while a combination of two fermionic or two bosonic 
quarks result in a bosonic field. Each sub-matrix in 
Eq. J7J) therefore consists of either fermionic of bosonic 
fields only. Although in principle arbitrary, the bosonic 
quarks are given the same masses as the corresponding 
valence quarks in order to cancel the contributions from 
closed valence quark loops. The above formalism is of- 
ten referred to as supersymmetric PQ^PT, although this 
only refers to the graded group structure of the matrices 
used in the construction of the theory. Furthermore, the 
term 'bosonic' merely indicates that those quark terms 
are treated as commuting variables. They are still spin 
1 /2 particles since the quark-antiquark pairs should build 
up particles with spin 0, representing the mesons of the 



theory. This violation of the spin-statistics theorem im- 
plies, as stated above, that in general the PQ theory is 
not a (fully causal) field theory. 

In the version of PQ%PT used in this paper, the sin- 
glet $o field has been dropped as discussed in Ref. 
This PQ singlet is expected to be heavy due to the axial 
anomaly in the same way as the if which is not included 
in unquenched xPT- Having defined the supersymmet- 
ric field matrix $ in Eq. (JJJ, one then proceeds by con- 
structing an effective Lagrangian for the PQ theory by 
requiring that this Lagrangian should be invariant under 
5t/(n va i + n sca |n va i) L x 5t/(nvai + n sca |n va i)_R in the same 
way as the Lagrangian of xPT was required to be invari- 
ant under SU (n /) l x SU (n /) R . This modification results 
in the same Lagrangian structure as for xPT, provided 
that the traces of matrix products in those Lagrangians 
are replaced by supertraces. The supertraces are defined 
in terms of ordinary traces by 



Str 



A B 
C D 



= Tr A - Tr D , 



(8) 



where A, B, C and D denote block matrices. For exam- 
ple, the block B corresponds to the [gygs] and [qsq~B] 
sectors of the field matrix in Eq. (JJJ and contains anti- 
commuting fields, while the block D represents the [gsgs] 
sector of Eq. J7J). Note that the removal of <&o implies the 
assumption 



Str ($) = . 



(9) 



The entire external field formalism introduced by Gasser 
and Leutwyler [2j for %PT can be generalized to include 
the extra degrees of freedom. In practice, external fields 
will be used in the valence sector only. 



B. Lagrangians and LECs 

We now proceed with the construction of the La- 
grangians for the Goldstone bosons and organize them 
according to the Weinberg power counting scheme. The 
external fields included are the vector and axial-vector 
fields and a^, as well as the scalar and pseudo-scalar 
external sources s and p. These are the suitably general- 
ized versions of those used in standard xPT 131 ■ Under a 
symmetry transformation g LtR e SU(n m i + ji scai \n vll i) LjR 
the fields transform as 



U 

X = 2B{s + ip) 
L 



9RUg[, 
9RX9 f L , 

9l l^gl-id^gr g[, 
9Rr tl g R -id fJ ,g R g R . (10) 

However, one can also define another set of quantities 
which behave differently under chiral symmetry trans- 
formations, and which turn out to be more useful for the 
construction of the Lagrangians. If one considers 



r u = v u 



- a,j 



exp 



(11) 
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then it is possible to find a matrix h such that 

u-> g R uh) = hug\. (12) 

One can then proceed by constructing a set of quantities 
that transform under chiral symmetry as X — > h X W . 
Such quantities are 

= i{v)(d ll -ir ll )u-u(d l _ L -il l _ L )v)} , 

zr = U Frv±^ u , (13) 

where Fx, and Fr denote the field strengths of the exter- 
nal fields I and r, such that F£ v = dH v - d u l» - i [P, i"]. 

is then defined analogously in terms of r. It should 
be noted that the Minkowski convention is used through- 
out this paper instead of the Euclidean one needed in lat- 
tice QCD, in order to maintain compatibility with the ex- 
isting literature in NNLO xPT. In terms of the quantities 
defined in Eq. I|13[l the lowest order, or 0(p 2 ) Lagrangian 
is given by 

C a = ?j-(u»u li +x+), (14) 

where the shorthand notation 

(A)=StrA (15) 

has been introduced. At this level there are two parame- 
ters F and B, which depend on the number of sea quark 
flavors in the PQ theory. In the power counting scheme, 
each derivative or factor of 1^ , r M counts as one, and each 
factor of s,p as two powers of the momentum p. The or- 
der in p is indicated by the subscript of the Lagrangian. 
The 0(p 4 ) Lagrangian has the generic form 

12 

£4 = LiXi + contact terms 

4=0 

= L Q {u^ifu^) + Li (u^Uf,) 2 + L 2 {u t "u u ){u t _ l u u ) 
+ L 3 (Kit M ) 2 ) + U (u fi u fJ ,)(x+) + U (u ll u ll x+) 

+ L 6 ( X+ ) 2 + L 7 ( X ^ + ^( X l+ X 2 _) 
+ iL u (x- (Vup-^C-U 

+ H^Fl+F^+^ixx^), (16) 

where the definition x- = X- — (x~)/ n sea has been ap- 
plied. Furthermore, the lowest order equation of motion 
is given by 

Vu P - \x- = 0. (17) 



The Lagrangian of Eq. i|16|) contains three types of terms. 
Of these, the terms proportional to Hi are contact terms 
which contain external fields only. Thus they are not rel- 
evant for low-energy phenomenology, but they are nec- 
essary for the computation of operator expectation val- 
ues. Their values are determined by the precise definition 
used for the QCD currents, and they are conventionally 
labeled h\ and H[ for unquenched xPT with rif = 2 and 
71 f = 3 quark flavors, respectively. The terms that de- 
pend on Ln and L12 are proportional to the equations 
of motion, and as such they can always be reabsorbed in 
higher order LEC:s, see Ref. f° r a full proof. The La- 
grangian at 0(p 6 ) is also known E3] and can be written 
in the form 

112 

Ce = KiYi + contact terms, (18) 

i=l 

where the form of the operators Yi can be found in 
Refs. 0,EJ. The PQxPT operators may again be ob- 
tained from the results of Refs. [2! fill if all traces are 
replaced by supertraces. All manipulations and identities 
used there to decrease the number of terms in the 0(p 6 ) 
Lagrangian, such as integration by parts and usage of the 
equation of motion, remain valid under this modification. 
This is why it is possible to apply the nj flavor results 
of Refs. 0, E 113 to PQxPT by using the appropriate 
graded matrices and supertraces. 

The divergences of PQxPT follow directly from those 
of rif flavor xPT calculated in Ref. [l6j, provided that 
rif is set equal to the number of sea-quark flavors. This 
is so since all the manipulations in Ref. E3 f° r 71 f na ~ 
vors remain valid when traces are replaced with super- 
traces. Alternatively, this derivation of the Lagrangians 
of PQxPT can be argued for with the Replica method of 
Ref. [l]| • The renormalization procedure is thus identical 
with the one of xPT also in the partially quenched case. 
An extensive discussion of this procedure can be found in 



TABLE I: The different sets of LEC:s for unquenched (n/) 
and partially quenched (n sea ) xPT- The number of physically 
relevant terms in each set is indicated by n p h, and the number 
of contact terms by n ct . The relationships between the various 
LEC:s are discussed in the text. 





XPT 
2 


XPT 
3 


X PT 

n 


PQXPT 
2 


PQxPT 
3 


LO 


F, B 


Fo,B 


jp(n) R (n) 


F, B 


Fo,B 


NLO 

n p h + n ct 


7 + 3 


LI 
10 + 2 


jr(n) 
11 + 2 


jj(1-pq) 

11 + 2 


11 + 2 


NNLO 

M p h + n c t 


53 + 4 


ci 

90 + 4 


112+3 


112 + 3 


112 + 3 
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Ref. 0| and references therein. The finite parts of the Lj 
and Ki for the different theories considered are summa- 
rized in Table J] For the case of xPT with two or three 
flavors the generic Lagrangians referred to above can be 
further simplified using the Cayley-Hamilton relations for 
2 x 2 or 3 x 3 matrices. These allow for a further reduc- 
tion of the number of operators, resulting in the different 
sets of constants shown in Table [I] The superscript r for 
the NLO and NNLO LEC:s indicates that they are the 
renormalized versions as defined in Ref. |16| . In particu- 
lar, the NNLO coefficients are given by the finite parts of 
the Ki, multiplied by F 2 in order to obtain dimensionless 
quantities. 

Consider next the LEC:s of PQxPT with n flavors of 
sea-quarks, which to lowest order are identical to those 
of unquenched n flavor %PT. At NLO, they are l^ npq ^ 
through L r ^2 Pq \ but as argued above the two terms 
jMjipq) an( j j^r(npq) ^ e removed us ing field redefi- 
nitions or the equations of motion 0, [Tg. For three 
flavors of sea-quarks, the standard choice is 



T r(3pq) _ T r(3pq) _ n 
^11 — ^12 — U : 



(19) 



whereas the unquenched two-flavor Lagrangian as defined 
in the first paper of Ref. 2| differs by an L\ Y type term. 
This makes no difference at NLO since that term does not 
contribute, but in order to get the correct correspondence 
at NNLO between the l\,c\ and the L^ (2m) , K- {2pq) one 
should take 



L 



r(2pq) 
11 



-11/4, L 



r(2pq) 
12 



0. 



(20) 



For the case of three sea-quarks, the L\ are simple linear 
combinations of the L^ 3p , which have been given in 
Refs. [ISini- Explicitly, 





_ jj(3pq) 


+ LS (3p9) /2, 




L r 2 


_ j^r(3pq) 


T r(3pq) 




LI 


T r(3pq) 
— ^3 


T r{3pq) 


(21) 



and 



_ r r (3pq) 

^4. ..12 — -^4. ..12 ) \ AL ) 

which indicates, as also discussed in Ref. that L^ 3p9 ^ 
is an independent LEC in partially quenched xPT with 
"-sea = 3, but in unquenched xPT it can be absorbed 
into L\, 17 2 and U z . One can therefore conclude that 



the numerical value of L 



r(3pq) 



cannot be determined by 



experiment, but is accessible via partially quenched Lat- 
tice QCD simulations. Similarly at NNLO, the C\ of 

unquenched xPT are linear combinations of the K^ 3 , 
and the corresponding relations have been derived in 
Ref. [13. 

The expressions given in the following sections of this 
paper concern PQxPT with n sca = 3 only, and therefore 



the superscripts (3pg) of the i^ 3p9 ^ and the K\ {apq> have 
been suppressed in most of the remaining equations. 



^r{3pq) 



C. Quark Masses and Propagators 

The version of PQxPT considered in this paper has 
three flavors of valence quarks (rival = 3), three flavors of 
sea quarks {n aca = 3) and consequently three flavors of 
bosonic 'ghost' quarks. For simplicity, the different quark 
masses are identified in the following calculations by the 
flavor indices i = 1,...,9, rather than by the indices 
u,d,s and V,S,B of Eqs. © and 0. The results are 
expressed in terms of the quark masses m q via the quan- 
tities Xi = ^Boiriqi such that X1jX2,X3j belong to the 
valence sector, \4, X5: X6 to the sea sector, and X7, X&> X9 
to the ghost sector. The latter ones do not appear in the 
results since the ghost quark masses are always set equal 
to the masses of the corresponding valence quarks, such 
that X7 = Xi,X& = X.2 and x.9 = X3- As discussed above, 
this is necessary in order to cancel the disconnected va- 
lence quark loops, which are replaced by the shifted mass 
contributions that involve the sea quarks. 

Furthermore, the quantities <i va i and d soa are used to 
indicate the number of nondegenerate quark masses in 
each sector. Thus <i va i = 1 has all valence quark masses 
degenerate, such that xi = X2 = X3j while d va i = 2 
indicates that xi = X2 7^ X3- Finally, in the case of 
dvai = 3, which is not needed in the present calculations, 
Xi 7^ X2 7^ X3- I n a similar way, for the sea quark masses 
4ea = 1 implies that X4 = X5 = Xe, while for 4ca = 2 
one has X4 — X5 7^ X6> an d for d sca = 3 all the sea 
quark masses are nondegenerate, such that X4 ^ X5 
Xe- The notation (e? va i + d sca _) is often used in this paper 
to indicate what degree of degeneracy in the quark masses 
is being considered. For example, the expression for the 
pseudoscalar meson mass for n va i = 2 and n SGa = 3 in the 
(1+2) case depends on one (degenerate) valence quark 
mass and two distinct (nondegenerate) sea quark masses. 

The meson propagators for the supersymmetric for- 
mulation of PQxPT have been calculated in Ref. and 
they correspond to the limit mo -> oo of the results in 
Ref. Q, where mo is the mass parameter of the singlet 
field $o- Again, for calculational reasons, the results of 
Ref. have been translated from the Euclidean formal- 
ism into Minkowski space. In general, three distinct types 
of propagators are encountered in the calculations. The 
simplest one is the propagator of a charged, or flavor-off- 
diagonal meson, connecting the meson field qiqj with its 
antifield qjc[i- It is given by Q 



iGUk) = 



k 2 ~ Xij + i£ 



(Mi)- (23) 



where the combination of quark masses Xij = (Xi + Xj)/2 
corresponds to the lowest order meson masses, and the 
signature vector ej is defined as +1 for the flavor indices 
of the rival + "sea fermionic quarks, and as —1 for the 
flavor indices of the n va i bosonic ghost quarks. In the 
present calculation, with the number of valence and sea 
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quarks as given above, ej assumes the values 

_ J +1 for j — 1, ... ,6 
6j ~ I -1 for j = 7,8,9. 



(24) 



The propagator of a neutral, or flavor-diagonal me- 
son, has a more complicated structure, since it connects 
mesons with different flavor indices as well. A propagator 
which connects the meson fields qiqi and qjq~j is written 
in the form 81 



(25) 



where the second term Gf^fc) is expressed in terms of 
a sum of single and double poles 0, • There are two 
distinct cases, as the double pole only appears if either 
i = 3 or Xi = Xj- For i ^ 3 and Xi^Xj, the propagator 
G\ d is given by 



iGUk) = 



R 



R 3 



k 2 



■ ie 



k 2 



D7T 



"J 



IE 



Xir + I? 



Xn 



- ie 



(26) 



where the factors R are referred to collectively as the 
single pole residues. For i = j or Xi = Xji the residues of 
the first two terms in eq. i|26|l become singular and the 
whole expression should be replaced by 



*G%{k) 



R'l 



{k 2 - Xi + ie) 2 k 2 - Xi 

_j_ TTlt 

k 2 - Xtt + i£ k 2 - x n + 



is 



- ie 



(27) 



where the residue of the double pole is denoted Rf. Note 
also the appearance of an auxiliary residue Rf in eq. (|27|l . 
All expressions for the propagator G?- depend on the 
lowest order neutral pion and eta meson masses in the 
sea quark sector, which are denoted by x-n and Xri- F° r 
d SC a, = 3 these are determined by the relations 



Xtt ~T Xt] 
XttXt] 



3 (X4 + X5 + Xe) , 

g 1 (X4X5 + X5X6 + X4Xe) , (28) 



which have no polynomial solution, but for d soa = 2 one 
has x-k = Xi an( i Xn = V^X4 + 2/3x6- The propagators 
for d sea = 2 can then be obtained from Eqs. (|2(j|) and l|27|) 
by taking the appropriate limits. In particular, the terms 
with a pole in x-n are 110 longer present. For d sea = 1, 
the fact that Xtt = Xij — Xi gives rise to several further 
simplifications, and in that case the terms with a pole in 
Xn vanish as well. 

The residues R of the neutral meson propagators in 
PQxPT are in general rational functions of the sea and 



valence quark masses which can be expressed in terms of 
the more general quantities R z a b defined by 



R-ab — 



^abc 



Xa 
Xa 



Xb, 
■ Xb 



R 



Xa ~ 
(Xa 



Xc 



Xb){Xa - Xc) 



abed 



R 



abedefg 



Xa Xd 
{Xa - Xb)(Xa - Xc)(Xa - Xd) 
(Xa ~ Xe)(Xa ~ X/)(Xa ~ Xg) ' 



(29) 



and so on. Note that R z a b has the same dimension as 
Xi for an even number of indices and is dimcnsionless for 
an odd number of indices. For the case of d soa = 3, the 
residues generated by Eqs. (|26|l and l|27(l are 



R 



3 hi 

Rf 

m 



Rii56jkl ) 
^2456^77 ) 



R 



Rl 



57T7] 



R 



Rl 



Rl 



(30) 



From these definitions, it is apparent that residues of the 
type Rf or Rj n „ vanish if i is a sea-quark index. Thus the 
propagators in the sea-quark sector of PQxPT contain no 
double poles as expected, since the origin of the double 
poles lies in the quenching of the valence quark loops. For 
the case of d sea = 2, due to the cancellations in the sea- 
quark sector (as discussed above) the residues simplify 
to 



Rjk 


— Rii6jk 


Rf 


— ^aion ' 


Ri 


= R\ n + 



R 



6)/ 



(31) 



so that the index tt no longer appears. For the case of 
^sca = 1) all residues associated with the sea quark sec- 
tor have reduced to numbers. Some nontrivial residues 
can still appear if the valence quarks are nondegenerate. 
These are 



j iij ' 



Ri = R 



ii I 



(32) 



where the double-pole residue has been retained mainly 
for notational consistency. Finally, it should be noted 
that if the sea quark masses are set equal to the valence 
quark masses, the propagator residues of PQxPT reduce 
so that the ir° and r\ meson propagators of unquenched 
xPT are recovered. 

Typically, a direct NNLO calculation with the propa- 
gators of Eqs. H26JI and (|27(l produces a large number of 
redundant residues in the output. This problem is espe- 
cially troublesome for the larger values of d sea and d V ai- 
It is thus useful to note that the various residues R sat- 
isfy a large number of algebraic relations, which provide 
an efficient, albeit tedious, way to simplify and compress 
the end results of the NNLO calculations. 
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D. Loop Integrals at NNLO 

The expressions for the NNLO masses and decay con- 
stants of the charged pseudoscalar mesons in PQxPT de- 
pend on several one- and two-loop integrals. After regu- 
larization and renormalization has been carried out, the 
finite contributions from these integrals are written in 
terms of the functions A, B and C, which are defined as 



Mx) = 



TTiexlogtx/V 2 ), 

Xt log(x»/M 2 ) - Xjlog(Xj/^ 2 ) 



-7T16 



Xi Xj 



-W(2x), 



(33) 



where fi denotes the renormalization scale and 7Tig = 
l/(16"7r 2 ). These integrals are often referred to as chi- 
ral logarithms, although the integral C does not contain 
any logarithm when all three arguments are equal. In 
the often encountered limit Xi — Xji the expression for 
B reduces to 



B(X,X;0) = -n 16 (1 + log( X /^ 2 )) 



(34) 



The functions described above are in principle suffi- 
cient to express all one-loop integrals encountered in the 
NNLO calculations, but the introduction of further com- 
binations of integrals is desirable in order to reduce the 
size and complexity of the results. For this purpose, the 
following three combinations of integrals have been in- 
troduced, 

A( X ;e) = A(x) 2 /(^i6X) 

+ ^ 16X (7r 2 /12 + l/2), 
B( X ,x;0,e) = A( X )B(x,x;0)/(n 16 x) 
- A(x) 2 /(2n 16X 2 ) 

+ ^l 6 (7T 2 /12 + l/2), 

B(xuXfAk) = XiB(Xi,Xf,V)+A(Xi), (35) 

of which the first two expressions are naturally generated 
by the dimensional regularization procedure, whereas the 
third one is useful since it is symmetric under the inter- 
change of Xi an d Xj- 

The NNLO calculation of pseudoscalar meson masses 
and decay constants also introduces a number of genuine, 
nonfactorizable two-loop integrals, which can be evalu- 
ated using a generalization of the methods described in 
Ref. |20|. The two-loop integrals encountered have the 
following general structure, 



(P0> 



d d q d d r 
(27r) d (27r) rf 



(36) 



X 



{q 2 - Xi) a {r 2 - X2) b ((q + r-p) 2 - Xd) c ' 



where a, 6, c = 1 or 2, and X represents the different com- 
binations of momentum factors q and r that can occur. 
The integrals thus generated by Eq. I|36|l are sometimes 



referred to as sunset integrals, and they can be expressed 
in terms of the H functions according to the definitions 
in Ref. [20| . The H functions satisfy a number of integral 
relations such that only H , Hi and if 21 are required to 
express all sunset integrals encountered. Furthermore, 
after regularization of the H functions, only the finite 
parts denoted H F , H F and H F i appear in the end re- 
sults. In the case of the decay constants, sunset integrals 
differentiated with respect to p 2 are also required, and are 
denoted by the primed quantities H F , H F and H F X . 

However, the appearance of double poles in the neutral 
meson propagators of Eqs. and l|27|l leads to a sig- 
nificant complication in Eq. (|36|l . as cases with a, 6, c 7^ 1 
will then show up. This added layer of complexity is 
accounted for by an extra integer argument n, which in- 
dicates the pole configuration of Eq. I|36|l according to 
Table [H] The finite parts of the H functions are thus 
generalized as 



H F (Xi,X],Xk;p 2 ) -> H F (n,Xt,X 3 ,Xk;p 2 ), 



(37) 



and similarly for the H F and H F \. In principle, eight 
different configurations can show up, and the value of 
n for each one of them is given in Table [HJ However, 
this expanded set of H functions obeys a generalization 
of the symmetries (under the exchange of mass argu- 
ments) discussed in Ref. |2(| , and thus some of the pole 
configurations turn out to be redundant, such that only 
n = 1,2,3,5 and 7 are required for a complete NNLO 
calculation of the pseudoscalar meson masses and decay 
constants. The case with n = 8 would however appear 
for a NNLO calculation of the coefficient of the double 
pole in the neutral meson propagator. Such a calculation 
has not yet been performed, but it should be noted that 
a NLO calculation of that quantity has been published 
in Ref. 0. 



TABLE II: Overview of the notation for the possible con- 
figurations of double poles in the H functions generated by 
Eq. 1361 in PQ X PT. Redundant configurations are given in 
parentheses. 





a b c 


n = 1 


1 1 1 


n = 2 
n — 3 
(n = 4) 


2 1 1 
12 1 
1 1 2 


n — 5 
(n = 6) 
n = 7 


2 2 1 
2 1 2 
12 2 


n = 8 


2 2 2 
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The explicit expressions for the expanded set of H 
functions with n ^ 1 can be obtained by differentiation 
with respect to the masses of the expressions for n = 1 
given in Ref. [2(J. It should also be noted that the finite 
contributions from the set of H functions contain inte- 
grals which have to be evaluated by numerical integra- 
tion, and these are again generalizations of the integrals 
for n = 1 in Ref. HJ. 

E. Notation for NNLO Results 

The length of the final NNLO expressions presents a 
problem. One major culprit has already been identified 
in the discussion of the PQ propagator residues, but even 
after that problem is dealt with, the size of the PQ ex- 
pressions far exceeds that of the analogous ones in un- 
quenched xPT. One reason for this is that the NNLO 
results in PQ^PT are highly symmetric under the in- 
terchange of quark masses, both in the sea and valence 
quark sectors, and thus contain a lot of unnecessary repe- 
tition. This observation suggests that the NNLO expres- 
sions could be efficiently compressed by summation over 
sea and valence quark indices. 

The summation conventions in the sea-quark sector are 
implemented through two new indices s and t, which can 
appear as indices of explicit quark masses Xi as well as 
among the indices of the residue functions R or the prede- 
fined combinations thereof. These sea-quark summation 
indices should be interpreted as follows: If an index t is 
present once or several times, then there will always be 
an occurrence of the s index as well, and the entire term 
is then to be summed over all pairs of different sea quark 
indices. If the index s is present but t is not, then the 
entire term is to be summed over all sea quark indices. 
Elementary examples are 

Xs = ^2 Xi, 

Xst = Xi i- ( 38 ) 

(y) = (45),(46),(56) 

Thus all the terms in the end result, where the depen- 
dence on the sea-quark masses is expressed in terms of 
the indices s and t, explicitly have the (required) sym- 
metry under interchange of sea-quark masses. It should 
be noted that the summation is over all three sea-quark 
flavors irrespective of the value of d sca . 

Further symmetries exist which can also be used to 
advantage in the compactification of the NNLO results. 
Firstly, the valence quark sector has a symmetry under 
the interchange of valence quark masses. This symmetry 
has been implemented by introduction of the summation 
indices p and q. In this paper, they are necessary only 
for the expressions with d va i = 2, and then always occur 
for the valence quark masses xi and X3- If the index q is 
present, there will always be an index p and the resulting 
sum is over the pairs (p,q) — (1, 3) and (p,q) — (3, 1). If 



only p is present, the sum is over the indices I and 3. As 
an example, consider 

Mx P )R p qn Xp = A(xi)R&,Xi + [1~3], (39) 

which demonstrates that any contribution written in 
terms of the (p, q) notation is symmetric (as required) 
under the interchange of the valence quark masses xi 
and X3- Secondly, for d sca = 3 the sea-quark sector 
exhibits an additional symmetry, under the interchange 
of the lowest order neutral meson masses x™ and Xv 
This symmetry has been implemented by the indices m 
and n. If the index m is present, there will always be 
an index n and the corresponding sum is over the pairs 
(m,n) = (tt, rf) and (m,n) = (t),tt). For example 

Mxm)KiiXm = A(xr,)KuXf, + [v <-> 4 (40) 

Similar to the earlier cases, if only the index n is present, 
then the term is to be summed over the x-n and Xn masses. 

The summation techniques described above already in- 
corporate most of the recurring combinations of terms in 
the end results, which has the added benefit of avoiding 
the introduction of an unreasonable amount of special- 
ized notation. However, certain combinations of quark 
masses x and propagator residues R appear through- 
out the results and have therefore been given dedicated 
names. As an example, combinations of the type 

Xg = 3 Yj Xi, 

i = 4,5,6 

Xt g = I E Ra biiXl (41) 

i = 4,5,6 

where the index g indicates the power of the sea-quark 
masses averaged, appear throughout the expressions. 
The corresponding ones for d sca = 2 can be obtained by 
setting — + R^ . Again, it should be kept in mind that 
the summation always runs over all three quark flavors, 
regardless of the value of d sea . The calculation of the 
pseudoscalar meson decay constant of Ref. ic| made use 
of many more complicated averages of the above type, 
with products of up to three residues. However, most 
of the quantities x which involve the residues R can be 
reexpressed in terms of simpler quantities, such that e.g. 

X7,a = - 1 , X7,i = "Xtt and xli = ~X V - The use of the 
X notation has therefore been discontinued, and the only 
quantities of the type given in Eq. H41[) that are used in 
this paper are Xi an< 4 Xi- The other named combina- 
tions consist of sums of products of quark masses and 
propagator residues. For c? sca = 3, these are 

t>v r>i i -pi o p« 

n ijkl — n jkk ' n jll ~ ZIt jkl , 

nu _ pi r>j pi pj 

^ijkl rt TTrik JX Tirik - n '7rr;fc- n "7r7;Z 

pi pj _i_ pi pj 

pum pi f? m — f?i F? m 

^ijkl ^Tnm'k m'jk ^ravri'k^m'jl 

pi pm _|_ pi pm 
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ryumn 



^m'ik^n' jk rt m'ik Jx n l jl 



R 



ijkl 



^m'il n'jl ' 

l '7l 



2R l v l R irikXkl + Kr,l R i v lXl, 



R 



ijkl 



rywmn 
-^ijkl 



r>i nm 
mm' k m' jkXk 



2R mm / faR m i jiXkl 



^^inm'l-^m' jkXkl T" ^rnm'l^m' jlXh 
^m'ik^ri jkXk ~r ^-H- m ' ik^ri jl Xkl 
^Rm'tlRn'jkXkl + Rm'ilRn'jlXh 



(42) 



and the analogous ones for d s 



2 are 



no 
^■ijkl 



Tjum 
-^ijkl 



Rl 4 + -R 



3 J 



kk 



2R 



jk > 



7? 



ijkl 



r>wm 
^ijkl 



^Vk^nk n r)k n -ql 

R ni RJ v k + R rii R nl' 

r>i rim r>i Dm 

^mk^jk ^mk^jl 

r>i T>m - m pm 
n ml n jk ' n ml lx jl > 

R \k R %kXk + 2R \k R \lXkl 

^Kjl^kXkl + H^l^iXU 

R mk Rl jkXk + 2 R mk R JlXkl 

2R l ml RJlxki + R m l R jlXl- 



(43) 



For the R u and i?™ type terms, the indices m, n denote 
either tt or 77 (independently of each other). If m denotes 
a tt then the corresponding index ml denotes an r] and 
vice versa. Similarly an n denoting a tt implies that the 
n' denotes an r\. As an example, consider 



^ijkl 



z>rj 
"■r/ik^TTjk 
pit nV 
^ml^TTjk ' 



' ^nik^wjl 



(44) 



which should suffice to explain the usage of the indices 



and 



Note that this particular combination is 



symmetric under the interchange of tt and ?y only if i = j, 
which means that this particular R u can appear in the 
final result only in combinations that are symmetric in 
the indices i and j. 

For the case of d sca = 2, the indices m', n' no longer ap- 
pear, since x-k = Xi an d has therefore canceled out from 
the expressions. In fact, the index m in Eq. (|43[) always 
represents an rj. However, Eq. (|43[l has been written in a 
slightly more formal way to illustrate the similarity with 
the expressions for c? SGa = 3. It is actually possible to ob- 
tain the expressions for d sea = 2 from the corresponding 
ones for d sca = 3 in Eq. Q42f) by removing all occurrences 
of the indices tt and ml there. It should be noted that 
the letters u,v,w in R U ,R V and R w are not indices, but 
rather indicate different types of residue combinations. 
The indices m, n used in the above formulas have noth- 
ing to do with the summation indices m, n defined earlier, 
which are encountered in the analytical expressions in the 
next sections. In Eqs. I|42|) and (|43|l they only indicate 
an occurrence of either a tt or an 77. 



The implementation of the flavor permutation symme- 
tries described in this section is a complicated task, as 
the direct output of a NNLO calculation of masses or de- 
cay constants using the propagators defined in Eqs. H26|) 
and (|27|l produces highly redundant expressions. Espe- 
cially for the higher values of d sea and d va i, the direct 
output sometimes consists of tens of thousands of terms. 
Thus the output has to be cleaned up using the many 
algebraic relations between the residues R, before any 
attempt at implementing the flavor permutation symme- 
tries can be made. This process involves the factorization 
of up to ~ 200 different expressions of varying length and 
complexity. As an example for d v& \ = 2, the parts of the 
direct output proportional to A(xn', e) and A(xi;e) alone 
are hundreds of terms long, while in factorized and sim- 
plified form, they contain only ~ 5 terms. After consid- 
erable trial and error in Maple, the simplification of the 
results was successfully carried through. In addition, the 
final simplified results have been checked algebraically 
against the original (long) ones. 



III. ANALYTICAL RESULTS FOR THE 
MASSES 

A. Masses at NNLO 

The corrections to the mass of a pseudoscalar meson 
are obtained by consideration of the self-energy contri- 
butions to the propagator of the interacting field theory. 
That propagator is defined in terms of the Fourier trans- 
form of the two-point Green's function, 



iA(p)= / dW^lT^^O^-p), 



(45) 



where $y = q^qj denotes any of the off-diagonal mesons 
in the valence sector of PQ%PT, and HI denotes the vac- 
uum of the interacting theory. To lowest order, iA(p) = 
G°j (p) of Eq. . When written in terms of the fields 
in the Lagrangian, the resulting self-energy diagrams can 
be summed as a geometric series [2l|. giving 



iA(p) 



Ml 



(46) 



where Mq denotes the lowest order mass of the meson 
which is being considered, and Xi m ^ denotes the de- 
pendence of the self-energy on all the lowest order me- 
son masses. The quantity Y,(p 2 ,Xi) takes into account 
the contributions from the one-particle-irreducible (1PI) 
diagrams. The physical masses, which include the inter- 
action, of the off-diagonal mesons in PQxPT are defined 
by the position of the pole in Eq. (|4fi|) . 



phys 



M 2 + £(M p 2 h v 4 ) 



(47) 



whereas the propagators of the neutral or diagonal 
mesons in PQ^PT have double poles even after resum- 
mation. 
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FIG. 1: Feynman diagrams up to 0(p 6 ) or two loops, for the self-energy E(Mp hys ,xO- Filled circles denote vertices of the £2 
Lagrangian, whereas open squares and diamonds denote vertices of the £4 and Ca Lagrangians, respectively. In the top row, 
the first diagram from the left is of 0(p 2 ), whereas the other two diagrams are of 0(p 4 ). The diagrams of 0{p 6 ), which give 
the NNLO correction to the meson mass, are shown in the bottom row. Of those, the third diagram from the left is called the 
"sunset" diagram in the text. 



The expression for the self-energy £ can be written 
as a string of terms which denote the 1PI diagrams of 
progressively higher order, 

S(M p 2 hys)X i) = S 4 (M p 2 hys)X4 ) 

+ E 6 (M p 2 hys , Xi ) + 0(p% (48) 

where £4 contains the self-energy diagrams of 0(p 4 ), and 
£6 those of 0(p 6 ). For the present NNLO calculation, 
the form of M 2 hyB has to be determined up to 0(p 6 ). It is 
sufficient, for this purpose, to use the lowest order mass 
instead of M 2 h in £ 6 since the diagrams in that term 
are already of 0(p 6 ). However, in the case of £4 the 
argument M^ hys has to be expanded in a Taylor series 

around Mq, since the diagrams in £4 are of 0{p A ). If 
M 2 hys is formally written as 



M, 



phys 



M£ + Mi 



Mi 



(49) 



then an expansion of £4 up to 0(p e ) gives 



£ 4 (M 2 



phy: 



, X i) = MK,Xi) + Ml 



dp 2 



(50) 



where M 2 represents the NLO correction to the lowest 
order meson mass, which is given by £4(^0, Xi)- The 
final formula for the pseudoscalar meson mass to 0(p 6 ) 
is thus 



^phys 



Ml + £ 4 (Af 2 , Xl ) 



£ 4 (M 2 , Xi ) 



dU(p\ Xl ) 



dp 2 

+ £ 6 (M 2 , X4 ) + 0(p 8 ), 



A A 



(51) 



where the last two terms are of 0(p 6 ) and represent the 
NNLO correction to the lowest order mass. The Feyn- 
man diagrams that contribute to ^e(M 2 ,Xi) are shown 
in Fig. In the case of the self energies at NNLO, the 
0(p e ) vertex introduces the LEC:s K"[ 7 through KT^, K^ 9 
and K r m . 

The physical mass of a pseudoscalar meson $jj is given 
to NNLO in the form 



M phys = Xij + 



£(4)vs 
1 



r(6)vs r(6)vs 

°Ct "loops 

F 4 



+ 0(p s ), (52) 



where the LO result Mq = Xij has already been inserted, 
and the 0(p 4 ) and 0(p 6 ) contributions separated. The 
NNLO contribution <5' 6 ) has been further split into the 
contributions from the chiral loops and from the 0(p e ) 
counterterms or LEC:s. The superscripts (v) and (s) in- 
dicate the values of d V ai and d sca , respectively. It should 
be noted that we have chosen to give the results to the 
various orders in terms of the lowest order decay constant 
Fq and in terms of the lowest order masses, since these 
are the fundamental inputs in PQxPT. The situation is 
different in standard x?T, where the main objective is 
comparison with experiment, in which case the formulas 
are most often rewritten in terms of the physical decay 
constants and masses. 



B. Expressions for d va i = 1 

In the previously published paper on the NNLO me- 
son masses in PQ%PT [9(, which treated only the most 
degenerate (1+1) mass case, an overall factor xi was fac- 
tored out from the analytical expressions. However, it 
has turned out that it is not possible to find such an 
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overall factor in a meaningful way for the less degenerate 
mass cases. This is especially evident for d va \ = 2. Thus, 
for consistency, the factor xi has not been factored out 
from the expression for the (1+1) mass case reproduced 
here either. That expression also uses the new R nota- 
tion as well as a more efficient set of loop integrals, but 
it is of course equivalent to the one presented previously 
in Ref. 0. 

The NLO result for c? va i = 1 is rather short, and there- 
fore it suffices to give this expression for d sea = 3 only. 
The corresponding results for d sca = 1,2 can easily be de- 
rived from that expression by consideration of the appro- 
priate limits, i.e. X5 -> Xi for d sca = 2 and X5 , Xs — *• X4 
for dgea = I- The combined NLO result (loops and coun- 
terterms), is 

,5(4)13 = _ 24 xi Xi ~SLlxl+ 48Z£xiXi 
+ 16 Llx\ - l/3A(xm)iCxXi 

- l/3A( X i)R c lX i 

- 1/3 B(xi, xi, 0)Rfxx, (53) 

which is in agreement with Refs. 0, Q . The NNLO re- 
sult has been split into the contribution from the loop 
diagrams and the contribution from the 0(p 6 ) countert- 
erms. The latter consists of the finite part of the right- 
most diagram in Fig. and is given for d soa = 3 by 



^ )13 = - 32 K[ rX l ~ 96K r 18X ixi ~ HUC&xl 

- 48X 2 r xiXi ~ 48^X2X1 

- 144^ 2 xi Xi ~ 16^X1 + 48^25X1 

+ K r 26 [96xiXi+48 X 2Xi] + 432^x1X1 
+ 32 K r 39 xf + 96/doXiXi- (54) 

As for the NLO result above, the contributions from 
the 0(p 6 ) counterterms for d sca = 1,2 can be derived 
straightforwardly by taking the appropriate mass limits 
of this expression. 

However, the situation is different for the loop con- 
tribution at NNLO, since it is much larger and has a 
rather complicated structure. This makes it a difficult 
task to derive the results for d sea = 1,2 directly from the 
rf soa = 3 case, and therefore the different cases are given 
separately below. The need for such expressions is clear 
since many Lattice QCD simulations work with d sca = 2 
rather than d sca . = 3. As expected, all infinities have 
canceled for all expressions in the renormalization pro- 
cedure, and the result for t^** is equivalent to the one 
published in Ref. Q. The FORM output with these ex- 
pressions can be downloaded from the website [14j, The 
chiral loop contributions to the pseudoscalar meson mass 
at NNLO in PQxPT, for d val = 1, are 



Cops = [3X1X4 + 26/3 XiX4-Xi] +4n 16 L r lX i + ""16 L r 2 [16 X iX 2 + 2 x?] + Trie £3 [3/2*1*4 

+ 17/3 XIX4- 5/2 Xl ] + n 2 16 [73/64 X1X4 + 15/32 X ?X4 - 3/32 Xl ] + 384i^XiX4 - 1152i^xiX4 

- 384L;LSx?X4 + 576 if X ixl ~ 384^x1X4 - 128 L^x? + 64ifxi - S A( X i) L r [ X \ 
+ R d lX i] +8l(xi)LIXi +20A(xi)^x? ~SA( X i)L r 3 [xl + Rfxi] + 16 A( X i) L\ X iXa 

+ A( X i)L r 5 [32/3 XT + 16/3 R d lX i] ~ Mxi) L r 6 [l6 X iX4 - 32 x 2 ] + 32 A( Xl ) L r 7 Rf X i 

- 64/3A(xi)^Xi + 5/9A(xi) 2 Xi + Mxi)B(xi, Xi, 0) [11/9*? + 1/9 JZ?*i] 

+ 2/9 A(xi)C(xi,Xi,Xi, 0)^X1 - A(xi,e)7r 16 [11/12 X \ ~ 1/4 12? Xi] + 3 A( Xu ) tt 16 X iX 4 

+ 24A(xi 4 )ioXiXi4 + 60 A(xi 4 ) £3X1X14 - 48 A{ X u) L r 5 X iXu + 96 A{ X u) L r s X iXu - 9/4 A( X i 4 ) 2 Xi 

- 2A(xi 4 )5(xi,Xi,0)xiX4 - A(xi4,e)vr 16 [9/2 X iX4 + 5/2 X \] + 128 A( X a) L\ X iXi 

+ 32 A{ Xi ) L r 2X iX4 - 128 Afa) £4X1X4 + 128 A( X a) L\ XiXa + 8/9 A( Xi )B( Xl , xi, 0) X iX4 

- 2A(x4,£)7n 6 xiX4 - 8B( Xl , Xl ,0)L r Rf X j - 8 B( Xl , Xl , 0) L\ Rf X \ + B( Xl , Xu 0) L r 4 [8 X ?X4 
+ 24i?fxiX4] + B( Xl , Xl ,0)L r 5 [8/3x? + 16i2?Xi] - B( Xl , Xl ,0) L r e [16x1X4 + 32i?fxiX 4 ] 

+ 16 5(xi,xi,0)^(^) 2 xi - B( Xl , Xl ,0)L r s [16/3 X \ + 32 R d lX \- 16/3 {Rff Xl } 

+ B(xi,Xi,0) 2 [2/9 Rixl + 1/18 [Rffxi] + 2/9 B( Xl , Xl , 0)C( Xl , X i, Xi, 0) {Rffxl 

+ 29/36 ^XiA^^x 2 + 16C(xi,Xi,Xi,0)£^i??x 2 X4 + 16/3 C( X u Xi, Xi, 0) L\ Rf X? 

- 32C(xi,Xi,Xi,0)^i?i ! X 2 X4 - 32/3C(xi,Xi,Xi,0)L^?X? + 5/9 H F (1, X i, Xl, Xi, Xi) x\ 
+ -ff F (l,Xi,Xi4,Xi4,Xi) [l/4xiX4~Xi] + 2i? F (l,xi4,Xi4,X4,Xi)xiX4 

+ 4/9# F (2,xi,Xi,Xi,Xi)#i < Xi + 3/4iJ i? (2,xi,Xi4,Xi4,Xi)^X 2 + 2/9 H F (5, X i, Xi, Xl, Xi) (Rffxl 

- 4 iff (3,Xi4,Xi,Xi4,Xi)-RiXi + 3/4^(1, xi,Xi4,Xi4,Xi)Xi + 6if^(l,X4,Xi4,Xi4,Xi)Xi 

- 3/4ff 2 F 1 (2,xi,Xi4,Xi4,Xi)i??X 2 , (55) 
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Cops = Tie £5 [-8/9 X^XiX4~X? + 26/3 xi X? + 35/9 X2X1] + 4tt 16 ^ X ? + ^6 L£ [22/3 Xl)X iX4 + 2 x? 

+ 26/3x2Xi] +^16^3 H/9 x^XiX4- 5/2 x? + 17/3 X1X1+ 43/18 X2X1] + ^ 6 [15/32 Xr ,XiX4 

- 3/32 x? + 15/32 xix?+ 43/64 x 2X i] + 384L^ XlX ? - 1152L^ x 2 Xl - 384L^Lg XlX 2 
+ 576L^ 2 x?Xi - 384 LgLg xi Xi - 128 L£Lg X ? + 64L£ 2 X ? - 8 A( Xr) ) L r Xr?Xl 

+ 16 A( Xv ) L r lX r,Xi + 4A(x„)£2Xr,Xi - 8 A(x„) £3^1X^X1 - 16A(x„)L; [Xr,Xi - Xi^ii Xi] 
+ 16/3 A( Xv )Ll [R v nX r,Xi+R V nXl] + ^A( Xv )L r 6 [x„Xi - Xi^i Xi] + 32 A( X „) ^ 46 i Xi 

- 64/3A(x,)^i?? 1 x 2 + 1/18 A(x v ) 2 m 2 Xi + 1/9 A( Xn )A( Xl ) R^R c lXl - 1/3 A( Xn )A( X46 ) HP n Xl 

- 4/9A( Xv )A( Xls )RlR z lsvXl + A( Xv )B( Xv , Xv ,0) [1/9x1^X1-5/18^x^X1] 

+ 2/9A(x,)S(xi,x t) ,0) [flJ^xi-flJXx*] + A( Xv )B( Xl , Xl ,0) [2/9 R^R^j + 1/9 RlXiXi 

- 1/27 Rl(Rl v ) 2 XiXs] + 2/9A( Xv )C( Xl , Xl , Xl ,0)R r > 1 R d lX 2 1 + A( Xv ,e)ir 16 [5/9 R^ X \ - l/4R c lXv Xi] 
-8A( X i)i5 [^X 2 + ^Xi] +8A(xi)^x? +20A(xi)^X 2 -8A( X i)^ [^X? + ^Xi] 

+ 16A(xi)^Xi^Xi + Mxi)L r 5 [32/3 R c lX 2 + 16/3 J2? X i] +^(xi)£@ [32 X ? - 16xi«f Xi] 
+ 32A(xi)L^fxi - 64/3A(xi)^i?Jx 2 + ^(xi) 2 [l/2 xi + 1/18 (i?J) 2 Xi] 
+ l/3A(xi)A(x46)i?l'iXi + 4/9A(xi)A(xi s )i?L^L ) Xi + 2/9A(xi)B(xi,X J ,,0)^ 1 ^X 2 
+ A(xi)B(xi,Xi,0) [ X \ + 2/9 (Rd 2 X 2 i + l/9^xi] + 2/9A(xi)C( X i,Xi,Xi,0)i?^x 2 

- A( Xl ,e)7r 16 [31/18 X 2 - 29/36^x 2 - l/4^xi] - ^(xi 4 ) 2 Xi - MxidMxm)xi - l/4A(xi 6 ) 2 Xi 
+ 48A( X4 )^XiX4 + 12A(x 4 )^XiX4 - 48 A( X4 ) L\ XlX4 + 48 A{ Xi ) L\ X iX4 

- 1/6 A{ X4 )B( Xv , Xti , 0)R v nX i X4 + 2/3A( X4 )B{ XuXv ,0)R n 14 Rl vXlX4 

+ l/3A(x 4 )S(xi,Xi,0)(i?l„) 2 XiX4 - 3/4A(x4,e)vri6XiX4 + 64 A( Xi6 ) L[ XlXi6 + 16 A( X 4s) L r 2 XiX46 

- 64A(x4e)i4XiX46 + 64A(x 4 6) ^6XiX46 + 4/9 A( X46 )B( Xv , Xv , 0) R^ XlX46 

- 4/27A(x4 6 )S(xi,Xr,,0) [RiUxl-RT&eXi] + 4/9 A( X4e )B( Xl , X i, 0) R^ XiX46 

- A(x46,e)7ri6XiX46 + -4(xis) Trie [1/4 XiXs +3/4xi Xi] + 8 A(xi s ) L r XlXu + 20 A( Xu ) L r 3 XlXu 

- 16A( Xls )L r 5XlXls + 32A( Xls )L r 8XlXls - A( Xls )B( Xl , Xv ,0) [4/9 RP U XlXls + 2/9 R v ls X 2 ] 

- 2/3A(xi s )5(xi,Xi,0) J R s 1 r; XiX s " A(xi s ,e)7r 16 [3/4 X iX, + 5/6 X 2 + 3/4 xi Xi] 

+ 8 B( XrnXri , 0)^x1^1X^X1 + 8/3 S(x„X„0) £5^4X1 - 16 ^(x^x^ 0)^X1^1X^X1 
+ 16B(x„X I) ,0)^(^ 461 ) 2 Xi - 16/3B(x„X„0)iS [#?i X 2 Xi - (^ 46 i) 2 Xi] 
+ 32B(xi,Xr ) ,0)i;^46i^Xi + 32/3 B( XuXv , 0) L r 8 R* 461 Rf Xl 

+ 2/9B( Xl , Xv ,0)B( Xl , Xl ,0)R 7 l 1 Rf X l - 8 B( Xl , Xl , 0) L r R* X \ - 8B( X i,Xi,0) L r 3 R d X l 
+ B(Xi,Xi,0)^ [8xi^X 2 +24xii?^Xi] + B(xi,xi,0) L£ [8/3 R\ X \ + 16 R d X 2 ] 

- B( XuXl ,0)L r 6 [16 Xl R c lX l+ 32 xi^Xi] + 16 B( X i, Xi, 0) L r 7 (R d ) 2 X i - B( Xl , Xl ,0) L r s [l6/3R c lX l 
+ 32i?^x 2 -16/3(^) 2 xi] +S(xi,Xi,0) 2 [2/9 R^Rf X \ + 1/18 (^) 2 Xi] 

+ 2/9B(xi,Xi,0)C(xi,Xi,Xi,0)( J Rf) 2 x 2 + 29/36 B{ Xu X i, 0, e) tt 16 i?f x ? 

+ 16C(xi,Xi,Xi,0)^Xi^X 2 + 16/3C(xi,Xi,Xi,0)L^fx? - 32 C( Xl , Xu X i, 0) L\ Xl Rf X \ 

- 32/3(7(xi,Xi,Xi,0)^^X? + 2/9^(1, x„X„Xi,Xi)(^'i) 2 X 2 + 4/9 ff F (l, x „ Xi, Xi, Xi) KiK xl 

- H F (l, Xr ,, Xu , Xls , Xl ) [1/3 R\ lX \ + 1/12 R^ u x v Xi] + H F (l, Xl , Xl , Xl , Xl ) [1/3 x? + 2/9 {Rl) 2 X 2 } 

- H F (l, Xl , Xls , Xls , Xl ) [1/2 R^xl- 1/4 ^X? + 1/12 Rfxi] + 3/4 ^(1, X i4, Xi4, X4, Xi) XiX4 

+ ff i? (l,Xi4,Xi6,X46,Xi)XiX46 + 4/9ff F (2,xi,X ?) ,Xi,Xi)^n^X 2 + 4/9 i/ F (2, X i, Xi, Xi, Xi) R\Rt Xi 
+ l/4 J ff F (2,xi,Xi s ,Xi„Xi)^X 2 + 2/9ff F (5,xi,Xi,Xi,Xi)(^) 2 X 2 

- 2/3fff(l, X r„Xi„Xi s ,Xi)^ s ^Lr,X 2 - 4/3 J fff(l,xi s ,Xi„Xi,XiX^Lr,X 2 

- 4/3fff(3,xi s ,Xi,Xi„Xi)^X 2 - l/4fl£(l,X„,Xi-,Xi.,Xi)^i a Xi 

+ H F (l, Xl , Xls , Xls , Xl ) [l/2i? s 1 „x 2 -l/4i?Jx 2 ] + 9/4i/ 2 F 1 (l,X4,Xi4,Xi4,Xi)x 2 
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+ 3i/ 2 F 1 (l, X 46,Xi4,Xi6,Xi)x? - 1/4/^(2, Xl , Xls , Xls , Xl )i^ X 2 , (56) 



Cop' = *™ L o H/9 X 7 rXr ) Xi-X? + 26/3 x lX ? + 35/9 X2 X i] + 4tt 16 ^ X ? + tt 16 L r 2 [22/3 X7rXr)Xl + 2 X? 

+ 26/3x2Xi] +7Ti 6 L5 [-8/9 x.X^Xi- 5/2 x? + 17/3 XiX? + 43/18 X2X1] + 4> [15/32 X7rX?)Xl 

- 3/32 x ? + 15/32 XiX?+ 43/64 X2X1] + 384L^ X iX? - H52L^x?Xi - 384L 4 Xg X i X 2 
+ 576 if x?Xi ~ 384 LJLS xi Xi - 128 L^ X ? + 64L£ 2 X ? - 8 A( Xm ) L r R™ u XmXi 

+ 16 A{ Xm ) L r lXmXl + 4A( Xm )L r 2Xm Xi ~ 8A( Xm )L r 3 Rl\ lXmXl - 16 A( Xm )L\ [ XmXl - X iKu Xi] 
+ 16/3i4( Xro )iS [^iiXmXi+^niiX?] + l6A( Xm )L r 6 [ XmXl - Xl R™ u X i] + 32 A( Xro ) L? i^ 456nl X i 

- 64/3i4(xm)i5iCnXi + l/18l(x m ) 2 « l n) 2 Xi + 1/9 A( Xm )A( Xl ) R^R c lXl 

- 4/9A( Xm )A( Xls )R" n \ s R z lsmXl + A( Xm )A( Xst ) [1/27 R% u R v mnst Xi + 2/27 R^ s R^ ls R z tssmn Xi] 

- 1/27 A( Xm )B( Xm , Xm ,0) (R™ ls ) 2 RZsXiXs - 2/27 A( Xm )B( Xm , Xn ,0) R^sRZsRlis XiXs 
+ A( Xm )B( Xm , Xl ,0) [2/9 (iCn) 2 X 2 - 2/27ii- ls i?™ s i? s 1 m „xiX s ] 

- l/27A( Xm )5( x „, x „,0)iCs(KuJ 2 Xi^ - A( Xro )B( X „, Xl) 0) [2/27iC.^Ci.^ m nXiX. 

- 2/9 i^ni&i X 2 ] + A( Xm )S( Xl , Xl ,0) [2/9K"ii^X 2 + l/9i? I T 11 ^Xi-l/27i?™ s (i 1 ' s 1 m J 2 XiX,] 
+ 2/9A( Xm )C( XllXl , Xl ,0)R™ 11 RUl + A( Xmi e)n 16 [-1/4 XmXl + 1/4 R™ u XmXl + 5/9 R% u X \] 
+ l/9A(x w )i4(x,)^ii^ n Xi - 8A( Xl )LS [^ X 2 + ^ X i] + 8A( Xl )^ x 2 

+ 20A( Xl )^ x 2 -8i4(xi)L5 [i??X 2 + ^Xi] + 16 A( X i) L r 4 Xl R c , Xl 

+ A( Xl )L r 5 [32/3 R c lX 2 1 + 16/3 R d lXl ] + A( Xl ) L r e [32 x 2 - 16 Xl ^ Xl ] + 32 A{ Xl ) L r 7 R d lXl 

- 64/3A( Xl )L^ x 2 + A( Xl ) 2 [1/2 Xl + 1/18 (i?J) 2 X i] + 4/9^( Xl )^( Xls )i?™ lsJ RL m Xi 

- l/27i4(xi)i4(x.t)(^,) 2iZ t«iiXi + 2/9A( Xl )B( Xm , Xl ,0)iCn^X 2 + A( Xl )5( Xl , Xl , 0) [ X 2 
+ 2/9(i?J) 2 X 2 + l/9i?^ X i] + 2/9A( Xl )C( Xl;Xl , Xl ,0)i?^ X 2 + A( Xl , e) tt 16 [-31/18 X 2 

+ 29/36 R c lX \ + 1/4 i?? Xl ] + A( X1s )tt 16 [l/4 XlXs + 3/4 XlXl ] + 8A( Xls )L£ XlXls 
+ 20A( Xls )L r 3XlXls - mA( Xls )L r 5XlXls + 32A( Xls )L r sXlXls - l/4A( Xls ) 2 Xl 

- l/2A( Xls )A( Xlt ) Xl - A( Xu )B( Xm , Xl ,0) [4/9i?™ lsXlXls + 2/9i?™ lsX 2 ] 

- 2/3A( Xls )B( Xl , Xl ,0)i? s 1 x „ X i Xs - A{ X u,e)n w [3/4 XlXs + 5/6 X 2 + 3/4 XlXl ] + 32 A( Xst ) L\ X iXst 
+ SA( Xst )L r 2XlXst - 32A( Xst )i4XiX,t + 32A( Xst )L r eXlXst + 2/9 A( Xst )B( Xm , Xm , 0) R^iKst XiXst 

- 2/27i4(x.t)B(Xm,Xi,0) [RlZt xl ~ R?Z Xi] + 2/27A( Xst )S( X7r , X?7 ,0) fi^?Xi 

- 2/27A( Xst ) J B( Xw , x „,0,fc)i?; is i?^ ls i? t z SSOT)Xl + 2/9A( Xst )S( Xl , Xl ,0)i? s 1 w i 1 , ^XiX,t 

- l/2A( Xst ,e)^ 16XlXst + 8B{ Xm , Xm ,0)L r 4Xl R™ uXmXl + 8/3 B( Xm , Xm ,0) L r 5 R™ llX 2 mXl 

- 16B( X 0)L5xi^ n XmXi + 16S( X rm X'irn 0) Z/7 (^m456nl) Xl 

- 16/3B( Xm ,x m ,0)i5 [«niiX™Xi-(^456„i) 2 Xi] + 32 B( Xm , xi, 0) L\ R z m456n iRf Xi 
+ 32/3S( Xm , Xl ,0)^^ 456nlJ Rf Xl + 2/9B( Xm , Xl ,0)B( Xl , Xl ,0)i?™ 11 i?f X 2 

+ 32B( X7r , Xj ,,0) J Lyi?^4 56 ^ 1 i?^ 4567rlX i 4- 32/3-B( XTr , Xr; ,0)L8i?^4 56 ^ 1 i?^4 56Tl X i - 8 B( X i, X i, 0) L T Q Rf X \ 

- 8B( Xl , Xl ,0)L r 3 R d x j + B( Xl , Xl ,0)L r 4 [8 Xl R\ X \ + 24 Xl R d Xl ] + B( Xl , Xl ,0) L r 5 [8/3 Rt X \ 
+ 16R d lX 2 } - B( XuXl ,0)L r 6 [16 Xl R c lX 2 1+ 32 Xl R d Xl ] + 16 B( Xl , Xl , 0) L r 7 (R d ) 2 Xl 

- B( Xl , Xl ,0)L r s [16/3 R{ X \ + 32 R d x \- 16/3 (i?f) 2 X ij + B( Xl , Xl , 0) 2 [2/9 fljitf x 2 + 1/18 (Rf) 2 Xl ] 
+ 2/9B( Xl , Xl ,0)C( Xl , Xl , Xl ,0)( J Rf) 2 X 2 + 29/36 B( XuXl ,0,e)ir 16 R d X 2 

+ 16C( Xl , Xl , Xl ,0)Ll Xl R d lX l + 16/3C( Xl , Xl , Xl ,0)L r 5 Rf X l - 32 C( Xl , Xl , Xl , 0) L\ Xl R d X \ 

- 32/3C( Xl , Xl , Xl ,0)L r 8 R d lX l + 2/9H F (l, Xm , Xm , Xl , Xl )(RZ 1 ) 2 xi 

+ 4/9ff F (l, Xm , Xl , Xl , Xl )iC n i^ X 2 - H F (l, Xm , Xls , Xls , Xl ) [l/3i?™ llX 2 + 1/12 R mn i s XmXl\ 
+ 4/9H F (l, X ^ XviXllXl )R; n Rl llX l + ff F (l,Xi,Xi,Xi,Xi) [l/3 X 2 + 2/9 (^) 2 X 2 ] 

- if F (l, Xl , Xls , Xls , Xl ) [1/2 i?L„X 2 - 1/4 Rtxl + 1/12 + l/2!Z F (l, X i s , X i t , Xst , X i) X i Xst 
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+ 4/9H F (2, Xl ,x mi Xi,Xi)KiiR d ixl + A/9H F (2, X i,Xi,Xi,Xi)R C iRfxi 
+ l/iH F (2,xi 1 Xis,Xis,Xi)R d 1 X 2 1 + 2/9 H F (5 :X i,Xi,Xi,Xi) (R d ) 2 xl 

- 2/3H[ (l,Xm,Xls,Xls,Xl) RnlsRlsmXl ~ 4 /3 iff (1 , \ls , Xls , Xl , Xl ) RnlsRlsm Xl 

- 4/3 H[ (3, Xls,Xl,Xls,Xl)Rf Xl ~ l/4i?^(l,Xm,Xls!Xls)Xl)-Rmnl S Xl 

+ ^(1, Xi, Xis, Xis, Xi) [1/2^X1-1/4^X1] + 3/2H F (l, X st,Xis,Xit,Xi)xi 

~ l/4i/ 2 F 1 (2, X i, X is,Xis,Xi)^X?. (57) 



C. Results for d va i = 2 

In general, the results for <i va i = 2 are similar in form 
to those for d va i = 1- However, they are about twice as 
long, since the number of independent combinations of 
quark masses that can appear in the propagator residues 
and loop integrals is much larger. As already noted for 
d V ai = 1, the expressions for d va i = 2 also become longer 
with increasing values of e? soa . The NLO expression for 
c? va i = 2 and d soa , = 3, which agrees with the result of 
Refs. 0,0, is 

^(4)23 = _ 24^xiXi3 ~ 8L£x?3 + 48igXiXi3 

+ 16^xi3 " l/3l(x P )i?^Xi3 

- l/3A(x ro )JC 13 Xi3, ( 58 ) 

while the contribution from the 0(p 6 ) counterterms for 
dva! = 2 and d sca = 3 is given by 

I 



C )23 = - 32iX[ 7 x? 3 - 96iq 8 XiX?3 " 8 tf^x^xia 

- 48iX 2 r xiXi3 - 48^X2X13 

- 144iX 2 r 2 XiXi3 - 16# 2 r 3 XiXi3X3 
+ 24#J 5 x*Xi3 + #26 [96xiXi3 

+ 48X2X13] + 432# 2 r 7 x?Xi3 + 32# 3 r 9 x? 3 
+ 96#I xiXi3- (59) 



Similarly to the case of d va \ — 1 , the expressions for other 
values of d sca can be straightforwardly obtained by taking 
the appropriate quark mass limits in the above expres- 
sions. Finally, the chiral loop contributions for d va i = 2 
at NNLO are, for d sea = 1, dsca = 2, and d sca = 3, re- 
spectively, 



<Cs = ^ L o [X1X13X3 + 3xi3X4 + 26/3 x? 3 X4- 2 x? 3 ] +4^ 6 ^Xi 3 +7r 16 L£ [16 XisX* + 2 X13] 

+ iT 16 L r 3 [7/6xiXi3X3 + 3/2xi3xl + 17/3xi3X4- 11/3 xi 3 ] + < [59/192 xiXisXs + 73/64 xisxl 

+ 15/32xi3X4 - 77/192 x? 3 ] + 384i^Xi 3 X4 - H52L^xi3xl - 384L^Xi 3 X4 + 5761^X13x1 

- 384L^Xi 3 X4 - 128L^x?3 + 64 x? 3 - MXv) *w [l/24x 9 Xi3 + 1/8 X13X4 - 1/4^X13X4 

- 1/12^X13] " Mx P )L r [4/3 Xp Xi3 + 16/3 ^x P Xi3+ 4/3^xi3] - A( Xp ) L r 3 [10/3x P Xi3 
+ 4/3 R p qXp Xi3 + 10/3^X13] + 16i(x P )^^Xi3X4 + MXp)LI [8/3 Xp Xi3 + 16/3 R p XpXl3 \ 

- 16A( x?) )^^Xi3X4 + A( Xp )L r 7 [16/3 R d pXpi + 32/3 R d pXl3 ] - A( Xp ) L r s [16/3 x P Xi3 + 32/9 R p X i X 3 
+ 64/9 R p q xl 3 + 8/9 (R d ) 2 } + A{ Xp f [5/36 X13 + 1/18 (f?£) 2 Xi3 - 1/54^ Xp ) + A( Xp )A( Xp4 ) [1/24 Xp 

- 1/8 x, + 1/24 RP x P + 19/72 R p Xq ] + A{ Xp )A( Xqi ) [1/4 X i 3 + 1/6 R p q Xq - 1/24 R« Xp + 1/72 R* Xq ) 
+ A( Xp )A( Xl3 ) [1/18 Xp- 1/36x13] - 4/27A(xp)A( X4 )x4 + A( Xp )B( XpyXp ,0) [11/54 f?£xpXi3 

- 1/27 RP q xl- 1/I8^xi3 + 1/54 R p R d p X P ] - Mx P )B(x q , Xq, 0) [l /54 R p q R d q Xq + 1/36 R d qXl3 ] 
+ 2/27A(xp)B(xi,X3,0)i?«XpXi3 + 1/27 A( Xp )B(xi, Xs, 0, k) R q p Xp + A{ Xp ,e)n l& [13/72 x P Xi3 
+ 5/72 x 9 Xi3 + 5/36 X i3X4 + 5/18 R p XpXl3 - 1/4 R p Xl3Xi ] - A{ XpA ) tt w [3/4 Xp Xi3 - 3/2 X i 3 X4 

- 3/4Xl3] + 12-4(Xp4)iSXp4Xl3 + 30A(Xp 4 )i5Xp4Xl3 - 24A(Xp4)^5Xp4Xl3 + 48 A(Xp 4 )^8Xp4Xl3 

+ 1/4 A(xp4)S(x 9 ,Xg, 0)^X13 - A( Xpi )S(xuX3,0) [4/9 x P Xi3 + 2/9 X i 3 X4 - 1/9 RP q X 2 p ] 

- A(xp 4 )5(xi,X3,0,fc) [2/9 xp 4 + 2/9x13- 1/9 R q pXq ] ~ A{ Xpil e)n lfi [9/4 X i 3 X4 + 5/4 Xl3 ] 

- A(xi)A(X3) [l/36xi3 + 1/54x4- 1/9 i?^?xi 3 ] + 8i(xi 3 )^Xi3 + 20 A( Xl3 ) L r 2 X 2 13 
+ 32i(xi 3 )^X?3 - l/4A(xi 3 ) 2 Xi3 + A(xi 3 )S(xi,X3,0) [l/9 XiXs + 4/9 X 2 3 ] 
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+ 4/9A( X i 3 )S(xi,X3,0,fc)xi3 - 31/18 A( X i3,e) ir 16 x 2 13 ~ 9/4 1( Xl4 )A( X 34) X13 + 128 A( X4 ) L\ Xl3Xi 

+ 32A( X4 )i2Xi3X4 - 128A(x 4 )^Xi3X4 + 128 A{ Xi ) L\ Xl3 X4 + 16/27 A( X4 )B( Xl , X 3,0)xisX4 

+ 8/27 A( X4 )B( Xl , X3 ,0, k) X4 - 2 A( X 4, e) tti6 Xi3X4 - B( Xp , Xp ,0)n 16 [1/24 R d p X9X13 + l/8i?pXi3X 4 ] 

- 4/3B(x P ,X P ,0)iS^pXpXi3 - 10/3B( Xp , Xp ,0)L r 3 R d XpXl3 + 8 B( Xp , Xp, 0) L\ RP q XpXl3Xi 
+ B( Xp , Xp ,0)L r 5 [8/3 BP qX l + 4/3^xiX3] - 16B(x P ,Xp,0)igi?Px P Xi3X4 

- B(Xp,Xp,0)£5 [16/3 ^x P + 8/3^X1X3] +B(Xp,Xp,0) 2 [1/8 X pXi3 - 1/72 R p q R d x, Xl3 ] 
+ 2/27B(x P ,Xp,0)B(xi,X3,0)i?^XpXi3 + 1/27 B( Xp , Xp , 0)S(xi, X3i 0, k) R^R d Xp 

+ B( Xp , Xp ,0,s)n 16 [1/4 R* Xp iXis + 5/36 R d pX \ 3 ] - 1/36 B( X i, Xi, 0)B( X3 , Xs, 0) RfRf X13 
+ 32/3 5(xi,X3, 0)^^X13 + 32/9 B(xi,X3, 0)^^xi3 + 16/3 B( Xl , X3 , 0, fc) L? R d R d 
+ 16/9B( X i, X3 ,0,k)L r 8 RfRi - H F (l, Xp , Xp , Xl3 , X i 3 ) [1/4 Xp Xi3 + 5/12 x 2 3 - 1/6 R p q X j 3 

+ 1/18 (i^) 2 X 2 3 ] + # F (l,Xp,Xi4,X34,Xi3) [1/4 x 2 3 + 3/4 J?«XpXi3 + 1/4 ^X«Xi3- 1/4^X13X4] 

- ff F (l,Xi,Xi3,X3,Xi3) [1/36 x 2 3 - 1/3 i?3^?X 2 3 ] + 1/4# F (1,X13,X13,X13,X13)X13 

+ 2ff F (l,xi4,X34,X4,Xi 3 )xi3X4 - ^ F (2,Xp,Xp,Xi3,Xi3) [l/18i?^x 2 3 - 5/36i^x 2 3 ] 

- ff F (2,x P ,Xi3,X 9 ,Xi 3 ) [1/18 R q p R d pX l 3 - 1/36 R d X \ 3 ] - 1/4 H F '(2 1 Xp; Xl4, X34, Xl 3 ) ^p Xp4Xl3 
+ 5/72if F (5,x P ,Xp,Xi3,Xi3)(i?p) 2 Xi3 + l/36!Z F (5,xi,X3,Xi3,Xi3)i??i?3Xi3 

+ fff(l,Xp,Xp,Xi3,Xi 3 ) [20/9 x 2 3 - 4/9 R p R q p x\ 3 ] - 2 H[(l, Xp , Xl4 , X34 , X13) R q P xl 3 

- 2ff 1 F (l,Xp4,X 9 4,X P ,Xi3)Xi3 + fff(l,Xi3,Xi,X3,Xi 3 ) [2/9 x 2 3 - 4/9 i^fl? Xis] 

+ 2/9^f(3,xi3,Xp,Xp,Xi3) [R p q R d pX 2 13 -R d pX 2 13 ] + 2/9 H[ (3, Xl3 , Xp , Xq , Xl3 ) R q R d pX 2 3 

- l/9ff F (7,xi3,Xp,Xp,Xi 3 )(i?p) 2 X 2 3 - 3/4ff F (l,Xp,Xp,Xi3,Xi3)x 2 3 

+ 3/4ff F (l,Xp,Xl4,X34,Xl 3 )-RpXl3 + 3/4^(1, Xp4,X«4,Xp,Xl3)^X 2 3 

- 3/4ff F (l,Xp4,X< Z 4,X 9 ,Xi3) [xl 3 'R q P x\ 3 ] + ^ 2 F i(l,Xi3,Xp,Xp,Xi3) [5/12 X \ 3 - 1/6 R p R q xl 3 ] 

- fl£(l,Xi 3 ,Xi,X3,Xi 3 ) [1/12 xi 3 - 1/3 i?3^Xi 3 ] + 3/4 X i3, Xis, Xis, Xis) X? 3 

+ 6ff F L (l,X4,Xl4,X34,Xl 3 )Xl3 - 3/4i7 F L (3,Xp4,Xp,X< Z 4,Xl3)-RpXl3 

- ff 2 F i(3,Xi3,Xp,Xp,Xi3) [1/6 R p q R d p X 2 i 3 - 1/12 R d p xU - #2 F i(3,Xi3,Xp,X 9 ,Xi 3 ) [l/6R q p R d p x\ 3 

- 1/12 R d x 2 13 ] + l/24ff F (7,xi3,Xp,Xp,Xi3)(^) 2 X 2 3 + 1/12 F F (7, X i3, Xi, X3, Xis) ^ X? 3; (60) 

<fis = ^16 £5 h8/9 X ,Xi3X4 + XiXi3X3- 2 x? 3 + 26/3xiXi3 + 35/9X2X13] + 4tt 16 ^x?3 

+ 7ri 6 £5 [22/3x„Xi3X4 + 2x?3 + 26/3x2Xi 3 ] + Trie £3 [-8/9x^X13X4 + 7/6 XiXisXs - H/3x? 3 
+ 17/3 X1X13 + 43/18 X2X13] + 7T 2 6 [15/32 Xr,Xi3X4 + 59/192 X1X13X3 - 77/192 x? 3 + 15/32 X1X13 
+ 43/64 X2X13] + 384L^xiX 2 3 - H52L^x 2 Xi3 - 384i^ X iX 2 3 + 576 LJ 2 xlxis 

- 384L^xiXi3 - 128£^X? 3 + 64L£ 2 X ? 3 + Mx P ) *w [1/4x1^X13-1/8x1^X13 
+ 1/12 R p qT) xl3- 1/24 R c p X q Xi3] - A( Xp )L r [16/3 R p qr] XpXl3 + 4/3 R c p XpXl3 + 4/3 R* Xl3 ] 

- A{ Xp )U 3 [4/3 J^XpXi3 + 10/3 R c pXp Xi3 + 10/3 R d Xl3 ] + 16^)^X1^X13 

+ A( Xp )Ll [16/3 R p qn x P Xi3 + 8/3 R c p x P Xi 3 ] - 16A( Xp )^Xi^Xi3 + A( Xp ) L r 7 [32/3 R d pXp 

+ 16/3 R d Xq - 8/3 (R d ) 2 ] - A( Xp ) LI [32/9 R p qv XlX3 + 64/9 R p qr) X \ 3 + 16/3 R c p XpXl3 + 8/9 (R d ) 2 ] 

+ A{ Xp f [1/8 X13 + 1/18 (i?^) 2 Xi3 + 1/54 R p qv R c p Xp - 1/54 {R c p ) 2 Xp + l/72(i?£) 2 xi 3 ] 

+ A( Xp )A( Xps ) [1/72 RP qn x P + 19/216 ^ Xq - 5/216 x P - 1/24 R% Xq + 1/27 R c p Xp ] 

+ A( Xp )A( Xqs ) [1/18 i?^ x 9 - 1/72 R%R z qsp xp + 1/216 ^^ sp Xq + 1/12 ^ Xl3 ] 

+ A(x P )A(Xn) [1/162 (^) 2 ^ s X, + 1/54 RP^ Xp - 1/90 R p qv R^ qq Xq + 13/270 R^ Xp 

+ 1/18 Xq - 1/135 x, - 1/54 R^R; xp + 1/216 R c p R v vl3 Xp + 1/72 R c p R^ 3 Xq 

- 1/36 R c p R v vl3Xri ] + A( Xp )A( Xl3 ) [1/18 xp- l/36xi3] - 1/18 A( Xp )A( Xi ) «,) 2 x 4 
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- A( Xp )A( X46 ) [2/27^X46- 1/9 R; p Xp] + A( Xp )B( Xp , Xp ,0) [11/54 R p qrj R c p XpXl3 - 1/27 R p qri R c pX 2 p 

- 1/18 R p qv R c pX 2 13 + 1/54 R p m R d Xp ] + A( Xp )B( Xp , Xv ,0) [2/27 R^R c p X \ + 1/27 ii^XiXs] 

- A(x p )B(x„X g ,0) [l/54i^i^x, + V36i^xi3] + 2/27A( Xp )B( Xl , X3 ,0)i^i^ XpXl3 

+ l/27^(x P )5( Xl , X3 ,0,/c)i?^ Xp + A( Xp ,£)tt 16 [-1/4X1^X13 + 1/8X1^X13 + 5/18 i?^x P Xi3 
+ 7/36 R c pXp Xi3 + 5/72^x g Xi3-l/72^Xi3] - ^(XpsKie [l/4 XpX i 3 - l/8 X i 3Xs - 1/4 Xia 

- 3/8xiXi3] + 4^( Xps ) J L5 Xj , sX i3 + 10A(x P s)L 3 x P sXi3 - 8^( Xps )L5 XpisX i3 + 16l( Xps )Lg XpsX i 3 
+ A( Xps )A( Xv ) [1/27 R; p Xp - 5/216 R; s Xp - 1/24 R^ Xq + 1/12 W m Xl3 - 5/72 R^ Xp 

- 19/216 RP qs R z psri x q + 1/72 i?? 3Xp ] - A( Xps )B( Xp , X?) ,0) [2/9 ^ XpXl3 + 2/27 iZJ, XpXs 
+ 1/27 ^ sX?Xs + 1/27 Rj 3 xl~ 1/27 i?? sX i X3 ] + 1/12 A( Xps )B( x „ x „ 0) ^ Xl3 

- A( Xps )B( Xl , X3 ,0) [1/27 R; vX 2 p + 4/27 Ri snXpsXl3 + 1/27 R% XlX3 ] 

+ A{ Xps )B( Xl , X3 ,0,k) [1/27 R q pr)Xq - 2/27 R% Xps - 2/27 R% Xl3 ] - A( Xps ,e)ir 16 [3/8 Xl3Xs 

+ 5/12 x 2 13 + 3/8 xixia] - l/2i( Xp4 )i( X96 ) X i 3 - A( x „)^ 16 [1/8 xi i^i 3 Xi3 ~ 1/24 R v vl3XvX i 3 

+ l/12R v vl3X 2 13 ] - A( Xn )L r [SR n 13XvXl3 + 4/3R v vl3XvXl3 ] + 1& A( Xr) ) L\ Xr)Xl3 + 4A{ Xr) ) L r 2Xr)Xl3 

- A( Xn )L r 3 [8i?? 3X?)Xl3 + 10/3i^ 13X?)Xl3 ] - 16 A( Xn )L r 4 [ Xj)Xl3 - Xl R\ 3 Xl3 ] 

+ A( X ,)iS [8/3 RP pp Xr ) Xis + 16/3 i?? 3 X 2 3 ] + 16^X^6 \XvX13 - Xi^ 3 Xis] 
+ A( Xn ) Uj [64/3 Rl XnXu - 64/3 Rj 3 X 2 13 - 4/3 (R* 4 ) 2 R v vl3 + 16/3 i?* 46p Xq ] 

- A(xr,) LI [64/9 i?? 3 XlX3 + 128/9 i?? 3 X \ 3 + 16/9 R" 13 X % + 8/3 i^ 13 XlX3 + 8/27 i^ 13 X ^ 

+ 16/27 R v vl3X 2 ] +A( Xn ) 2 [1/54 Xl ^ 13 + 1/18 (i?? 3 ) 2 Xl3 - 1/54 i?^ 13X , + 1/54 i?^ 13Xl3 

- 5/108 R v vl3 Xv - 1/54 (^ 13 ) 2 X , + 1/72 (i^ 13 ) 2 X i 3 ] - 1/36 A{ Xv )A( X4 ) R v nl3 Xi 

+ A( Xv )A( X46 ) [1/9 Rj 3Xl3 + 2/27 R v vl3X46 ] + A( Xv )B( Xp , Xp ,0) [l/9 R^R^ X 2 p ~ 1/27 R; p R d pXp 

- 1/54 RP qq R d p Xp + 1/216 R d p R» nl3 Xp + 1/72 R d R^ l3 Xq - 1/36 R d R^ l3 Xv \ 

- A( Xv )B( Xp , Xv , 0) [2/27 R p m R^ XpXv + 1/27 R^R^ XqXv - 2/27 R p m R\ 3 X 2 p - 1/27 R^R^ XlX3 ] 
+ A{ Xr ,)B{ Xr ,, Xr ,,0) [1/9 X ii?? 3X i3- 5/18 i2? 3 x,Xis] + 2/27A( X7) )S( Xl , X3 ,0)i?^i?^ Xj)X i3 

+ 1/27 A(x v )B(xi,X3,0, k) R^R^Xv + A{Xr,,e)^w [-1/4 x^Xia + 1/8 Xi^isXis + 1/4^x^X13 
+ 5/9^ 3X 2 3 + l/8i?- 1 3 X , Xl 3 + 5/36i?- 1 3 X 2 3 ] - A( Xl )A(xs) [l/54R p v R qXq ~ 1/9 R\ v R% X i3 
+ 1/36 i^Xis] + 8A( Xl3 )^ x 2 3 + 20A( Xl3 )^ x 2 3 + 32 A( Xl3 ) L r 6 x 2 13 - l/4i( X i 3 ) 2 X i 3 
+ A( Xl3 )S( Xl , X3 ,0) [l/9xiX 3 + 4/9 xls] + 4/9A( Xl3 )B( XuX3 ,0,k) Xl3 - 31/18 A( Xl3 , e) tt 16 X 2 3 

- i( X i 4 )A( X34 ) X i 3 - l/4A( Xl6 )A( X36 ) Xl3 + 48i( X4 )L^ X i 3X 4 + 12A( X4 )L^ X i 3X 4 

- 48A( X4 )L r 4Xl3X4 + 48A( X4 )L r eXl3X4 + A( X4 )B( Xp , Xv ,0) [2/9i??X4X P X4 + l/9i?^ 4X9X4 ] 

- 1/6 A( X 4)B{xr l ,Xr l , 0)^3X13X4 + 2/9 A{ X4 )B( X i , X3 , 0) Rl v R 3 4r] Xl3X4 

+ 1/9 A( X4 )B{xi,X3,0,k) Rl^R^Xi - 3/4i( X4 ,e)7r 16Xl3X4 + 64 A( X46 ) L\ Xl3X46 
+ 16^( X46 )L2Xi 3 X46 - 64 A( X46 ) L 4X i 3X46 + 64A( X46 )Lg X i 3X46 

+ ^( X46 )S( Xp , Xv , 0) [2/81 R u p ; 46 XpXn - 2/81 JQj, XpX ^ - 2/81 i?^ 46 XlX3 + 4/81 iQ, Xp 
+ 2/81 i?;; 46 X J - 2/81 A(x46)B(x P , X„ 0, k) [R^ i6 Xp + R% 46 Xp ] 

+ 4/9A( X46 )B( Xv , Xr] ,0)R r > 3Xl3X46 - A( X46 )B( Xl , X3 ,0) [2/81 iZ^xiXs - 4/81 R? 346 Xl3 ] 

- i( X46 )B( Xl , X3 ,0,fc) [4/81^5*346X13 - 2/81^348] - A( X46 , £) 7T 16 X1SX46 

- B(x p , Xp ,0)tt 16 [l/8 Xl i?^ Xl3 + l/24i?^ XgXl3 ] - 4/3B( Xp , Xp ,0)L5i?^ XpXl3 

- 10/3B( Xp , Xp ,0)L r 3 R d pXpXl3 + 8B{ Xp , Xp ,0)L r 4Xl R p qvXpXl3 

+ B{ Xp , Xp ,Q)Ll [16/3 i?^ X 2 X i 3 - 8/3 i?^ X i X i 3X3 ] - 16 B{ Xp , Xp , 0) L r 6 Xl R p n XpXl3 

- B(x P ,Xp,0)L r 8 [32/3 i?^ X 2 Xl3 - 16/3 i^ X i Xl3X3 ] +^( Xp , Xp ,0) 2 [l/8i?^ XpXl3 

- l/72^^ X?Xl3 ] + B( Xp , Xp ,0)B( Xp , X)) ,0) [2/27i?? 3 ^ X 2 + l/27i?^ XlX3 ] 
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+ 2/27B( Xp ,Xp,0)B(xi,X3,0)R q pr ,Rix P Xi3 + 1/27 B(x P ,X P ,0)B(xi,X3,0,fc) R q pr ,R d p X P 
+ B( Xp ,Xp, 0,e)7r 16 [l/8 Xl R d Xl3 + 7/36R d X P Xi3 + 5/72 R d XqX i 3 ] 

+ B( Xp ,x v ,0)L r 7 [32/3 R d p R z i6qXp + 16/3 R d p R z i6qXq ] + B( Xp , Xn ,0) L r 8 [32/9 R d R z 46q Xp 
+ 16/9 R d R z i6qXq ] + 8B(x„X„0)L;xii2? 3 X,Xi3 + 8/3S(x J) ,Xr ) ,0)^^ 3 x'xi3 

- 16 B(x^,x^ 0)^X1^3X^X13 ~ 8 B( Xv ,Xv, 0) L r 7 RP 13 (^ 4 ) 2 Xi3 - B(x,,X,.0)iS [16/9 i?? 3 X13X4 
+ 32/9i?^ 3 xi3xi] - l/36B(xi,Xi,0)B(x 3 ,X3,0)i?^Xi3 + 32/3 B( Xl , Xs, 0) RfR d 3 Xis 

+ 32/9B(xi,X3,0)L5i2?i2g X i3 + 16/3 B( X i, X3, 0, fc) L? jtfiig + 16/9S(xi,X3,0,fc)Lg^ 

- H F (1, Xp , Xp , Xi3, Xis) [1/4 XpXis + 3/8 X ? 3 + 1/18 {R%?xlz - 1/6 R p qv R L p xls + 1/24 (i?pf X? 3 ] 
+ H F {1, x P , xi., Xss, X13) [-5/12 i?^ x 9 Xi3 + 1/6 R p qv XpsX i 3 + 1/8 x 9 Xis + 7/24 R p v Xl3Xs 

- 1/12 R c pXpsXl3 ] + l/2i/ F (l,x P 4,X 9 6,X46,Xi3)xi3X46 + H F (1, Xn , Xp, Xi3, Xis) [2/9i?^ 3 x 2 3 

- 1/18 i?^^ 13 x 2 3 + 1/9 RP pp R c p xl 3 + 1/36 i^isXis] +^ F (l,X„X^Xi3,Xis) [2/9 K 3 ) 2 X? 3 
+ 2/9 i?? 3 i?- 13 x 2 3 + 5/72 (i?^ 13 ) 2 x? 3 ] - tf F (l,X„,Xis,X3s,Xis) [1/8 R? ps XpXis + 5/24 R% s XqXl3 

- 1/3 i?? 3 x?3 + 1/24 R v vpsXnX i3 + 1/24 ^ 13 Xi3Xs] - # F (l,Xi,Xi3,X3,Xi3) [l/18i?^x 2 3 

- 1/3 Rl n R 3 lvX 2 13 - 1/36 R1R C 3X 2 13 ] + l/4ff F (l,xi3,Xis,Xi3,Xis)x 2 3 

+ 3/4ff F (l,xi4,X34,X4,Xi 3 )xi3X4 - # F (2,Xp,Xp,Xi3,Xis) [l/18i?^x 2 3 - 5/36i?^x 2 3 ] 

+ H F (2, Xp , Xri , Xl3 , Xl3 ) [1/9 R v pp R d p xl 3 + 1/36 R d p R v nl3 xi 3 ] - H F (2, Xp , Xl3 , Xq , Xl3 ) [l / 48 R%R d X \ 3 

- l/36R c q R d pX 2 3 ] - l/12H F (2, XpiXls , X3s , Xl3 )R d XpsXl3 + 5/72 H F (5, Xp , Xp , Xl3 , Xis) (i?p) 2 X? 3 
+ l/36i/ F (5,xi,X3,Xi3,Xi3)i?fi?3X 2 3 + H F (l, Xp , Xp , Xl3 , Xl3 ) [2x 2 3 + 4/9(i?^) 2 x? 3 

- 4/9i?^^x 2 3 + 2/9(i?p 2 x 2 3 ] + 2/3Hf(l,x P ,Xi s ,X3 S ,Xi3)^^ P X 2 3 

- 2/3 iff (1, Xps , Xqs , Xp , X i 3 ) R p sri X 2 i3 - 2 /3 H i( l ,X P s,Xqs,Xn'Xi3)Ri 3 R Z s PV X 2 i3 

+ 2/9iff(l,Xi3,X P ,Xr ) ,Xi3) [R^R^xls-RlpRqpnRpXls] - H F (l, Xl3 , Xri , Xri , Xl3 ) [2/9i?? 3 ^ 13 x 2 3 
+ 1/9 (^i 3 ) 2 X 2 3 ] + #f(l,Xi3,Xi,X3,Xi3) [2/9i?^x 2 3 -4/9^i?? )) x 2 3] 

+ 2/9i/ 1 F (3,xi3,Xp,Xp,xi3) [R p qv R d P xl 3 -R c P R d P xl 3 ] + 2/9H[(3 lXl3 , XpiXq , Xl3 )R; v R d pX 2 13 
+ 2/9H F ('i, Xl3 , Xp , Xv , Xl3 )Rl 3 R z pq71 R d pX 2 l3 - l/9i/ 1 F (7,xi3,Xp,Xp,Xi3)(^) 2 X 2 3 

- 3/4^(1 

, Xp? Xp, Xi3> X1 3 )X 2 3 - 1/4^1(1 , Xp, Xis, X3s, Xl3) R q r,R Z sqp Xl3 

+ i/4H F (i, Xps , Xqs , Xp , Xl3 ) [R p mX \ 3 + R% X \ 3 - R c pX \ 3 ] + l/4H F (l,Xp S ,X 9S ,X 9 ,Xi3)i?^i?: M X 2 3 

+ l/4-ff F l(l,Xps,X9s,X»?,Xl3)-Rl3-Rp W -Rsp);Xl3 ~ I/ 4 #21 (1 , Xrj , Xis , X3s , X13) R^Rll-qRld-n Xl3 

+ ff F (l, xis, Xp, Xp, Xis) [3/8 x 2 3 + 1/6 (R p qv ) 2 X 2 13 - 1/6 R p qv R c p X \ 3 + 1/24 (^) 2 x? 3 ] 

- ^(^Xis^X^Xis) [1/6 ^i?- 13 x 2 3 - 1/12 ^ 13 X 2 3 ] + l/24^ 2 F 1 (l,Xi3,Xr J ,Xr ) ,Xi3)(^i 3 ) 2 X?3 

- # 2 F i(l,Xi3,Xi,X3,Xi3) [1/6 R p qv R c qX 2 13 - 1/3 ^i??,X 2 3- 1/12 i??i?3X 2 3 ] 

+ 3/4^(1, Xl3,Xl3,Xl3,Xl3)Xl3 + 9/4 J ff F 1 (l,X4,Xl4,X34,Xl3)Xl3 + 3 / 2 #2l(l, X46, Xp4, Xg6, X13) Xl3 

- l/4i/ 2 F 1 (3,Xps,Xp,X 9 s,Xi3)i?pX 2 3 - #2 F i(3,Xi3,Xp,Xp,Xi3) [l/Q R p qri R d X \ 3 ~ 1/12 R c p R d pX 2 13 ] 

- H F (3, Xl3 , Xp , Xq , Xl3 ) [1/6 R; v R d pX 2 13 ~ 1/12 R c q R d pX 2 13 ] + 1/12 H F (3, X i 3 , X P , X„ Xis) i?^ 13 X? 3 

+ l/24i/ F (7,xi3,Xp,Xp,Xi3)(^) 2 X 2 3 + l/12# 2 F i(7,Xi3,Xi,X3,Xi3)i?^X 2 3, (61) 



Cops = Tie ^ H/9 X.X.X13+X1X13X3- 2 x? 3 + 26/3 xiX 2 3 + 35/9x2X13] +4tt 16 ^x?3 

+ n 16 L r 2 [22/3x^X13 + 2 x? 3 + 26/3x2Xi 3 ] + tt 16 [-8/9x^X13 + 7/6 XiXisXs - 11/3 Xis 
+ 17/3 xiX?3 + 43/18 X2X13] + 7T 2 6 [15/32 XxXr,Xi3 + 59/192 XiXisXs - 77/192 x? 3 + 15/32 XlX \ 3 
+ 43/64 x2X13] + 384L^xiX 2 3 - H52L^x 2 Xi3 - 384L^x lX 2 3 + 576if XiXis 
- 384L^xiX 2 3 ~ 128L^X? 3 + 64 X ? 3 + ^(Xp)^i6 [1/4x1^X13-1/8X1^X13 
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+ 1/12 EP^xXz- 1/24 # pXgX i 3 ] - A(Xp)L r [16/3 ^XpXis + 4/3 i^XpXia +4/3 i^xis] 

- A( Xp )L r 3 [4/3^x^3 + 10/3^x^X13 + 10/3^X13] + 16 A( Xp ) L£ Xl i^ Xl3 

+ A(x P )L r 5 [16/3 i?^ XpXl3 + 8/3 ii pXpXl3 ] - 16 A( Xp ) L£ Xl i^ Xl3 + A( Xp ) i r 7 [32/3 i^ Xp 

+ 16/3 R d Xq - 8/3 (i?^) 2 ] - A(x P ) LI [32/9 XlX3 + 64/9 RP qnrj X 2 13 + 16/3 R c p XpXl3 + 8/9 (R d p ) 2 ] 

+ A( Xp ) 2 [1/8 xi3 + 1/18 (i?^) 2 Xi3 + 1/54 R p qnri R c p Xp - 1/54 (R c p ) 2 Xp + 1/72 (^) 2 X i 3 ] 

+ A( Xp )A( Xps ) [1/72 RP q „ v Xv + 19/216 R^ Xg - 5/216 Xp - 1/24 Xq + 1/27 fl pXp ] 

+ A(x P )A( Xqs ) [1/18 R^x, - 1/72 RP vri R z qsp x P + 1/216 RP S7ir] R z qsp Xq + l/12i$xi 3 ] 

+ ^(Xp)^(Xm) [-l/27-R^ m „-R™g Xp + l/QRq mn Rnl3Xp + ^ l^l RP qmn R™ 13 X q + 1/27 R P qmn R mn \3 Xp 

- 1/54 R p qmn R v mnl3 Xm + 1/162 (i?? m „) 2 i?™ s X s ~ 1/54 i?™,^ X P + 1/216 R c p R v mnl3 X P 
+ 1/72 R c p R v mnl3Xq - 1/36 i^iCmaXm] + A( Xp ),4( Xl3 ) [1/18 Xp - 1/36 Xl3 ] 

+ A( Xp ),4( Xst ) [1/162 R u ppst x P - 1/81 i^ st Xp + 1/81 X P - 1/81 R^t X P - 1/162 R^ pst ] 
+ A( Xp )B(xp,Xp,0) [n/54RP q ^R c pXp Xi3 - 1/27 R^ n R c p X 2 - 1/18 i?^^ X 2 3 + 1/54 B^B* Xp ] 
+ A( Xp )B(x P ,x m ,0) [2/27 R™ 13 R c pX 2 p + 1/27 R™ 13 R c p XiX 3 ] ~ A( Xp )B( Xq , Xq , 0) [1/54^^, 
+ l/36i^ Xl3 ] + 2/27 A( Xp )B(xuX3,0) R 9 p« n RpXpXu + 1/27 A( Xp )B( X i,X3,0,k) R q pnr ,R c pXp 
+ A(x P , e) tt 16 [-1/4 Xl i?^ xia + 1/8 Xi^p Xis + 5/18 i?^ XpX i 3 + 7/36 R c p XpX i 3 + 5/72 R c p Xq Xi3 
-l/72^ Xl3 ] +A( Xps )^ 16 [-l/4 XpX i 3 + l/8 X i 3Xs + l/4 x 2 3 + 3/8 X iXi 3 ] + 4^( Xps ) L r oXps Xi3 
+ 10A(xps)L 3 x P sXi3 ~ 8A(x ps )LlxpsXi3 + 16 A{x ps ) LI XpsX i 3 - 1/4 A( Xps )A( Xgt ) Xi3 
+ A(x ps )A(xm) [1/27 R™ p Xp - 5/216 i?™ sXp - 1/24 iZ^X, + l/12i?™ gXl3 - h/72R™ qsXp 

- 19/216 R™ 8 R z psmXq + 1/72 iCi 3 X P ] + A( Xps )B( Xp , Xm , 0) [-2/9 i?™ s XpXl3 - 2/27 R™ qsXp Xs 

- 1/27 i?™ s x 9 Xs- 1/27 iC 13 x P + 1/27 KisXiXa] + 1/12 i( Xps )B( x „ Xg , 0) i? ? d Xl3 

- A( Xps )B( Xl , X3 ,0) [1/27 i?^ X 2 + 4/27 i?^ XpsXl3 + 1/27 i?^XiXs] 

+ A( Xps )B(xi,X3,0,k) [1/27 R^ n x q - 2/27 Ri nri xps- 2/27 R^xu] - A( Xps ,s)tt 16 [3/8 Xl3Xs 
+ 5/12 x 2 i3 + 3/8 xiXis] + ^( X m)^i6 [ -l/8 X ii?^ nl3X i 3 + l/24i?^„ 13XmX i 3 - l/12R v mnl3 x\ 3 \ 

- A{xm)L r [8-R™ 13XmX i 3 + 4/3i?^„ 13XmX i 3 ] + 16^4( Xm )L^ XmX i 3 + 4A(x 

- A(xm)L r 3 [8i?™ 13XmX i 3 + 10/3 iC„i 3X mXi 3 ] - 16A(x m )L r 4 [XmXl3 - X l-?Cl3 Xl 3 ] 

+ A( Xro )L£ [16/3 iC 13 XmXi3 + 16/3 i?™i 3 X? 3 + 8/3 iCmaXmXia] + 16A( Xm )Lg [ XmXl3 

- Xi^Cis Xis] + A( Xm ) L r 7 [32 R™ 13 XmXi3 - 32/3 i?™ 13 XlX3 - 64/3 i?™ 13 X 2 3 - 16/3 R d R z mi56np 

+ 16/3 RqR Z mi5&np - 8/3 (-R^456np) 2 + 16/3 R Z m imn\R Z mi5&n3\ + MXrn) L r s [-128/9 i?™13 XlX3 

- 64/9 i?™ 13 X 2 3 - 16/9 RpRmittnp + 16/9 -Rg-Rm4 56 „ p - 8/9 (i?m4 56np ) 2 + 16/9 i?m456„l-Rm456n 3 

- 8/3 Cnl3Xm - 8/3-R^„ 13X i X3 ] + A{Xm) 2 [1/81 RP mn Rn ps R™ ss Xs - 1/&1 Rl mn Rnp S RnssXs 

+ 1/18 (R^ 3 ) 2 xi3-1/54RZ 3 R mnl3 Xm 

+ 1/54R™ 3 R 

Xi 3 + 1/162 (R™ ss ) R v mn i 3 Xs 

- 1/54 (R v mnl3 ) 2 Xm + 1/72 (R v mnl3 ) 2 Xl3 ] - A{ Xm )A{xst) [l/27 R™ 13 R v mnst Xl3 + 1/81 R^ st Xq 
+ 1/162 R v mn i 3 R v mnst Xm + 1/81 R p ™ t - 1/81 R p ™ t + 1/162 R p ™™ - 1/81 Ri 3 ™™] 

+ A( Xm )B( Xp ,Xp,0) [1/9 R p qmn R™ pp X 2 p - l/MR™ p R d pXp - 1/27 R™ 3 R d pX p - 1/72 R d R v mnl3 x P 
+ 1/72 R d R v mnl3Xq - 1/36 R d p R v mnl3 Xm] + A( Xm )B( X p,X m ,0) [-2/81 R p smn R™ S R™ S s x P Xs 

- 1/81 R P mn Rn qs Rnss X q Xs + 2/27 RnppRnl3 X P + V 2 7 i?™ p i?™ 13 X 1 X3 ] 

+ A(xm)B(x P , Xn, 0) [-2/81 R p smn RZ s R mqs XpXs - 1/81 R p smn RZ s R mqs X,Xs 

+ 2/27 R™ pp R n ml3 x 2 p + 1/27 iC pp iC 13 XiXs] - 1/27 A( Xro )B( Xm , Xm , 0) i?™ ls i?™ 3s i?™ s Xl3Xs 

- A( Xm )B(x mi Xn-; 0) [l/81 R n p S R nss R"mps XpXs H~ 2/81 R n p S R nss Rmqs Xl3Xs\ 
+ 1/81 A(xm)B(x m?Xn?0,fc) [R nps R nss R mps Xs R n psRnssRm q s Xs] 

- 1/27 A( Xm )B(Xn,Xn,0) RZ s R m lsR m 3sXl3Xs + A( X m )B( X l , X 3 , 0) [2/27 R q pmn Rl l p p X pXl3 
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- 2/81iC s i? s 1 m „^L„Xi3X s ] + A( Xm )B(xi,X3,0,k) [1/27 R q pmn R™ p X P ~ 1/81 ^.RLAn X.] 

+ A{Xm,s)lT W [-l/4XmX!3 + 1/8 XlR V mnl3 Xl3 + 1/4 Kl3 Xm X 13 + 5/9 i?™ 13 X? 3 + 1/8 R mn 13 X>nXl3 

+ 5/36iC„isXi 3 ] + MXn)A( X r,) [1/270 XiK,i3^i3 + 1/45 i?™ 13 K m i 3 Xn + 1/270 i?™ 13 K m i 3 Xi3 

- 1/81 R nss R r ip S R^ r p S Xs + 1/81 R nss R r ip S R^ r g S X s 29/1080 R^ppR^pp Xp + 1/72 R^ pp R!^ pp Xq 

- 1 7/540 R; pp R* qq xi3 + 1/5 ^13^13X13- 1/162 R* ss Rl ss R v mnl3 Xs] - A( X i)A( X3 ) R^ n R c qXq 

- 1/9^^X13 + 1/36^X13] + 8A(xi 3 )^Xi3 + 20 A( X is) X13 + 32A( Xl3 )L£ x 2 3 

- l/4i(xi 3 ) 2 Xi3 + A(xi 3 )B(xi,X3,0) [1/9X1X3 + 4/9 X ? 3 ] + 4/9 4(xis)B(Xi, X3, 0, k) X13 

- 31/18 i(xi3,e)7ri 6 x? 3 - l/4A( X i s )A(x 3s )xi3 + 32 A( Xst ) L\ X i 3Xst + 8 A{ Xs t) L r 2 X i 3 Xst 

- 32A(x st )L r 4 xi3Xst + 32A(x st )Lgxi3X,* + A( Xst )B(x P , Xm, 0) [1/81 X P Xm ~ 1/81 XpX™ 

- 1/81 R^ XiX3 + 2/81 R p ™ t x P + 1/81 R pq ^t Xq] - l/8lA( X st)B( Xp ,X m ,0,k) [R™ tXp 

+ R pq n st X P ] + 2/9A( Xst )B(xm,Xm,0)i?™ 13 i?™ t xi3X st + A{x s t)B{x„x^) [1/81 iQ X**> 

- 1/81 ^ g 7 t X.Xr ) + 1/81 iC3x P + 2/81 iC?Xis] +% i )%„X„0,fc) [-1/81 i^Xp 

- 2/81 R^ t X13 - 1/81 ^ + 1/81 - i( Xst )B(xi,X3, 0) [l/Sl i^ 3st XlX3 - 2/81 i?r 3st Xis] 

- A(x rt )B(xi,X3,0,A:) [2/81 i^ 3st Xl3 - 1/81 i?^*] - 1/2 e) tt 16 X i 3Xst 

- B(x p ,Xp,0)tt 16 [l/8xii?pXi3 + l/24i?^x 9 Xi 3 ] - 4/3 B( Xp , X P , 0) L r R d Xp Xi3 

- W/3B( Xp ,x P ,0)L r 3 R d p x P Xi3 + 8B( Xp , Xp ,0)^ Xl i^ XpXl3 - B( Xp ,x P ,0) L r 5 [16/3 XpXis 

- 8i?^4 Xl3 ] - 16B( Xp , Xp ,0)Lg X i^XpXi3 + B( X p, X p,0)L^ [32/3 R p q7zr] XpX 2 3 - 16 RP q7VT] X 2 X i 3 ] 
+ B(xp,Xp,0) 2 [l/8i?^XxpXi3-l/72i?^^ X9Xl3 ] + 5( Xp , Xp ,0)B( Xp , Xm ,0) [2/27 iCi 3 # P Xp 

+ l/27iCi 3 *£xiX3] + 2/27B( Xp , Xp ,0)B( Xl , X3 ,0)i?^i?^ XpXl3 

+ l/27S(xp,Xp,0)S(xi,X3,0,A)i^ r ,i^Xp + S( Xp , Xp , 0, e) tt 16 [l/8 Xl i?^ Xl3 + l/8i^ XpXl3 

+ 5/36i?p X 2 3 ] + £?( Xp , Xm ,0)L7 [32/3 i? p i?^ 456ng Xp + 16/3 iZ p iZ m456 „ gXg ] 

+ #(Xp,Xm, 0)-Lg [32/9 i?pi?^ 456n(3 Xp + 16/9 i?pi?^ 456n(jX(? ] + 8 B( Xm , Xm, 0) L\ X i-R™i3 X mXi3 

+ 8/3B( Xm , Xm ,0)^i?™ 13X 2 reX i 3 - 16B( XrojXro ,0)L£ X iiCi 3 XmXi3 

+ 16-B( X m, X m,0)L7i?^ 456nl i?^ 456 „ 3X i 3 - -B( X m, Xm , 0) Lg [l6/3 i?™ 13 XmX i 3 

- 16/3 -Rm 456 „i-R m456 „ 3 Xl3] + B(Xtt,Xv^) L 7 [I 6 / 3 -R^456r;p^4567rpXp + 32/3 i?^. 456j)p i?^ 4567r9X i 3 ] 

+ -B(X7r,X)j7 0)Lg [l6/9i?^ 456j?p i?^ 456ffp Xp + 32/9i?^ 456j)p i?^ 456x q X i 3 ] 

- 16/3 B(xtt, Xn-> 0' k) Uj [Rl 456vp R^ 4567rp — Rl 456rip R* 4567Tq ] — 16/9 B{x-k, Xr) , 0, fc) Lg [i?* 456j)p i?* 4567rp 

- i?; 456w ^45 6 J - l/36B( Xl , Xl ,0)B( X3 , X3 ,0)i?fi? 3 d Xl3 + 32/3S( Xl , X3 ,0)L r 7 i?^ Xl3 
+ 32/9B( Xl , X3 ,0)^i?fi? 3 d Xl3 + 16/3B( Xl , X3 ,0,fc)^i?fi? 3 d + 16/9 B( X i, Xs, 0, k) L r s RfR d 

+ H F (l, Xp , X P , Xis, Xis) [-l/4XpXi3 - 3/8 X 2 3 - 1/18 {R p ^fx\ 3 + 1/6 R p ^R c p X ? 3 - 1/24 (R c p fx 2 13 ] 
+ H F (l, x P , Xis,X3s,Xi3) [-5/12 X?X i 3 + 1/6 R p qnr) Xp sXi3 + 1/8 R p nr , x q Xi3 + 7/24 R p nr) Xl3Xs 

- 1/12 i? pXpsX i 3 ] + l/4# F (l, Xps , X9 t, Xst , X i 3 ) X i 3XiS t + # F (l, Xm , Xp , X i 3 , X i 3 ) [2/9 R p mn R™ 13 X 2 i3 

- 1/18 RP qmn 

ran 13 Xl 3 + 

1/9 R n pp 

Rp X13 + 1/36 -R p 

''ranlS x? 3 ] 

+ H F (1, Xm , Xm, Xi3, Xis) [2/9 (iCi 3 ) 2 X 2 3 + 2/9 iC 13 iC„ 13 X? 3 + 5/72 (R v mnl3 ) 2 X 2 3 ] 

+ ff F (l, xm, Xis, X3s, Xis) [-1/8 R™ s x P Xi3 - 5/24 i?" s Xq Xi3 + 1/3 i?™ 13 X ? 3 - 1/24 R v mnps XmXl3 

- l/24iC„i 3 Xi3X s ] + ^ F (l,X.,X„Xi3,Xi3) [l/9i?; pp i?^ pX 2 3 + 2/9i?^ 13X 2 3 

+ l/36iC,i3^Ci3X 2 3 ] + ^ F (l,Xi,Xi3,X3,Xi3) [-l/18i?^ X 2 3 + l/3i?LAx 2 3 

+ l/36i?^ X 2 3 ] + l/4if F (l, Xl3 , Xl3 , Xl3 , Xl3 ) X 2 3 - H F (2, Xp , Xp , Xl3 , Xl3 ) [1/18 R p q7rv R d p x 2 13 

- 5/36R C pR d xl 3 ] + H F (2, X p,Xm,Xi3,Xi3) [1/9 B%„R* X \z + 1/36 R d p R v mnl3 xis] 

- H F (2, X p,Xi3,X q ,Xi3) [1/18^,^X13-1/36^X13] - 1/12 H F (2, Xp , X i„ X3s , X i 3 ) i?^ XpsX i 3 
+ 5/72i/ F (5, Xp , Xp , Xl3 , X i 3 )(i?^) 2 X 2 3 + l/36i/ F (5, Xl , X3 , Xl3 , X i 3 )i?fi? 3 d X 2 3 
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+ iff (1, Xp , Xp , Xi3, Xis) [2 x? 3 + 4/9 (^) 2 X? 3 " 4/9 R p q ^R c p x? 3 + 2/9 (R c p ) 2 X 2 13 \ 

+ 2/3H[(l > Xp,XU,X3 S ,Xl3)Rqnri RZ sq P Xl3 ~ 2/3 H l 0-> Xps, Xqs, Xp, X13) R^ v *13 

- 2/3 fff(l, Xps, Xgs,Xm,Xl3)-R™13-RspmXl3 + 2 /9 fff (1, Xl3, Xp, Xm, X13) [^mn-RmnU Xl3 
R npp R qpm R p X13] ~ H?(l,xi3,Xm,Xm,Xia) [2/9 ^nl3^nml3 Xl3 + 1 /9(Kinl3) X13] 

- 2/QH[ (1,X13,X7T,X77'X13) [K?,ppRq n pK,nl3 Xl3 ~ ^13^13 X13] 

+ ^f(l,Xi3,Xi,X3,Xi3) [2/9i?^x?3-4/9^i?? m7 x? 3 ] + 2/9 iff (3 

! Xl3; Xp; Xp; Xl3) [-^q-Tr^-^p 

X?3 

- R p R i X13] + 2/9 iff (3,Xl3,Xp)X«>Xl3)-R^r n -R p Xl3 + 2/9 iff (3, Xl3,Xp,Xm,Xl3) R7l3 R pqm R t Xl3 

- 1/9 iff (7,xi3,X P ,Xp,Xi3) (Rp) 2 X 2 i 3 - 3/4iff(l,Xp,Xp,Xi3,Xi3)x? 3 

- l/4if2 F i(l,Xp I Xi s ,X3 S ,Xi3)-R^^pX?3 + l/4/f 2 F i(l,Xp Sl X9 S ,Xp,Xi3) [-Rg^Xis + -R?^Xi3 - -RpXis] 

+ l/^ifaiC^XpsiXgsjXg^ls) Xl3 + 1 / 4H 2l( 1 1 Xps 1 Xqs-; Xm-> X13) ^nlZ^pqm^ spm Xl3 

- l/4if 2 F i(l,Xm,Xls,X3„Xl3)i?™l3^1 m ^s3mX?3 + (1 , Xl3 , Xp, Xp , Xl3 ) [3/8*?3 + 1/6 (i^,) 2 *^ 

- 1/6 f?^x 2 3 + 1/24 (i?^) 2 X 2 3 ] - #2 F i(l,Xi3,Xp,X m ,Xi3) [1/6 i?? mn K nnl3 X 2 3- 1/12 i?^„ 13 X 2 3 ] 

+ l/24iff(l,Xl3,X>n,X>n,Xl3) (-Rmnl3) 2 Xl3 + V I 2 #21 (1 ) Xl3 , Xtt , Xn , Xl3 ) -R^13^irl3 Xl3 

+ ^ 2 F i(l, X13, Xi, X3, Xis) [-1/6 R p q „ n R c q x 2 3 + 1/3 Rl^RU x 2 3 + 1/12 i2fii§ x? 3 ] 

+ 3/4iff(l,Xl3,Xl3,Xl3,Xl3)Xl3 + 3/4^(1, Xst, Xps, X?t,Xl3)Xl3 

- l/4iff 1 (3,Xp,,Xp,X 9S ,Xi3)^X 2 3 - # 2 F i(3,Xi3,Xp,Xp,Xi3) [1/6/^/^X13 - 1/12 i^Xis] 

- -fff(3,Xl3,Xp,Xg,Xl3) [1/6 -Rp^-Rp Xl3 - 1/12 -Rg-RpXl3] + l/12-ff F l(3,Xl3,Xp,Xm,Xl3)-Rp-Rmnl3Xl3 

+ l^fff^T.xia.Xp.Xp.XiaJ^Xia + V 12 Hftf, X13, Xi, X3, X13) RfR d 3 Xia- (62) 



IV. ANALYTICAL RESULTS FOR THE DECAY 
CONSTANTS 

In the previous section, explicit analytical results for 
the masses of charged pseudoscalar mesons to NNLO in 
PQxPT were presented for all possible degrees of de- 
generacy in the quark masses. The decay constants of 
the pseudoscalar mesons are similarly known, and have 
for the most part been published earlier in Ref. l n 
that reference, analytical results were presented for all 
cases except for the most difficult one with d va \ = 2 and 
dsca = 3. However, the more advanced simplification 
methods used for the meson masses in this paper has 
also made it possible to bring that expression down to 
a manageable size. This is not surprising, since the re- 
sults for the decay constants are in general slightly less 
complex than the ones for the meson masses. In order to 
obtain a complete set of NNLO calculations, the expres- 
sion for the decay constant with d va i = 2 and d sca = 3 is 
presented here. 

The decay constants F a of the pseudoscalar mesons are 
obtained from the definition 



(OK(O)|0 o (p)> = ivV*f«, 



(63) 



in terms of the axial current operator A%(0). In the fol- 
lowing developments, the flavor index a has been sup- 
pressed for simplicity. The Feynman diagrams that con- 
tribute to the axial current operator at NNLO, or 0(p 6 ), 
are those shown in Fig. [21 Diagrams of 0(p 2 ) and 0(p 4 ) 



also contribute to Eq. (|63|) via the renormalization of the 
pseudoscalar meson wave function <fi a (p)- From Eq. 163fl . 
the expression for the decay constant of a pseudoscalar 
meson to 0(p 6 ) is given by 



F, 



phys 



Z 



= Fn 



F 4 (M 2 



phys 



Xi) 



F 6 (M p \ ys ,Xi) + 0{p% 



(64) 



where V Z is the wave function renormalization factor. 
Here the subscripts of the matrix elements F indicate 
the chiral order, and should not be confused with the 
flavor index a in Eq. i|63|) . Thus F4 and Fq denote the 
matrix elements of the axial current operator at NLO and 
NNLO, respectively. In the above equation, the lowest 
order contribution F2 has already been identified with 
Fq. The wave function renormalization is given in terms 
of the self-energy diagrams by 



z- 1 = 1- 



5S(p 2 ,x.) 



dp 2 



(65) 



M- 



which becomes, when expanded such that all contribu- 
tions up to 0(p 6 ) are taken into account, 



E' 

Z = 1 + T 

dE 4 (p 2 ,x,) 



-E' 2 



dp 2 



d£ 6 (p 2 ,x.) 



M- 



dp 2 



(66) 
(67) 



M- 
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FIG. 2: Feynman diagrams up to 0(p r ') or two loops, for the matrix element F(Mp hyB , \i) of the axial current operator A^(0). 
Filled circles denote vertices of the £2 Lagrangian, whereas open squares and diamonds denote vertices of the £4 and C& 
Lagrangians, respectively. In the top row, the first diagram from the left is of 0(p 2 ), whereas the other two diagrams are of 
0(p 4 ). The diagrams of C(p 6 ), which give the NNLO correction to the decay constant, are shown in the bottom row. 



and is thus seen to contain self-energy contributions of 
0(p 4 ) as well as 0(p 6 ). It should be noted that the two 
terms in the expression for £' are of chiral order p 2 and 
p 4 , respectively. That the self-energy of 0(p 8 ) does not 
contribute at NNLO is evident from Eq. 1)64(1 . since the 
matrix clement of the axial current operator itself is of 
0(p 2 ) or higher. By combining all the above results and 
expressing everything, analogously to Eq. (|5 II) . in terms 
of the lowest order meson mass M , the final result for 
the decay constant at NNLO is obtained as 



C(p 4 ) contribution 



-Fphys = Fq + Fi(Xi) + Fr 



2 dp 2 



(68) 



Mi 



+ Fr 



3/ 9S 4 (p 2 ,Xi 



+ F 4 ( Xi 



dp 2 



Fr 



2 dp 2 



Mi 



Mi 



Mi 



2 dp 2 

F 6 (M 2 , X i) + 0(p s ) 



where all noncontributing terms have been discarded. In 
principle, the NLO result indicated by braces in the above 
equation depends on M 2 hys and has to be expanded, anal- 
ogously to the calculation of the meson masses, around 
M 2 up to 0(p 6 ). However, the extra terms generated by 
such an expansion can be eliminated using the identities 



dF 4 {p 2 , X t 
dp 2 



0, 



d 2 £ 4 (p 2 ,X«) 
d{p 2 ) 2 



0. 



(69) 



The former one follows from the fact that the matrix 
element F4 does not depend on p 2 , and the latter one is 
valid since the highest power of the momentum p that 
appears in the self-energy £4 is p 2 . The NNLO decay 
constants of the charged pseudoscalar mesons so obtained 
depend on the 0(p 6 ) LEC:s K{ 9 through K% 3 . 



The physical decay constant of a charged pseudoscalar 
meson $ij is given to NNLO in the form 



-Fphys — Fq 



f(4)v! 

~~F 2 ~ 



,(6> 
Jet 



+ // 6) 



loops 



F 4 



o( P 8 



(70) 

where the 0(p 4 ) and 0(p 6 ) contributions have been sepa- 
rated and denoted by / rather than 8 in order to minimize 
the potential for confusion with the rather similar ex- 
pressions for the pseudoscalar meson masses. Again, the 
NNLO contribution has been further split into the 
contributions from the chiral loops and from the 0(p 6 ) 
counterterms. As for the meson masses, the superscripts 
(v) and (s) indicate the values of d va i and d sca , respec- 
tively. Up to NLO, the pseudoscalar meson decay con- 
stant for (ivai = 2 and d soa = 3 is given by 



(4)23 



\2L\xi + 4££xi3 

+ A( Xp ) [1/6 BP^- 1/12 R;] 

+ l/4A( Xps ) - 1/12 A( Xm )RZ 

- l/l2B(x P ,x P ,0)K 



mn 13 



(71) 



and agrees as expected with Refs. 0- The LEC:s at 
0(p 6 ) give the following contribution to the pseudoscalar 
meson decay constant for d va i = 2 and d sca = 3, 



f (6)23 _ 4K r 2 

Jet — 4 1X 19 Xp 



2AK r 2oX iX.i3 + 24^X2 



72i^ 2 x 2 + 8^ 3 xiX3. 



(72) 



The expression for the chiral loop contribution is rather 
long, but it can be brought to a manageable form with 
the new notation introduced earlier in this paper. The 



result for d 



val 



2 and eL 



3 is 
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/loops = Tie L r [4/9 XvXr, - 1/2 XiX 3 + xi 3 - 13/3 XiXis - 35/18 X2 ] - 2 tt 16 LJ x? 3 - *16 ^2 [H/3 X.X„ + X? 3 
+ 13/3x 2 ] + ^16^3 [4/9 x.X^ - 7/12xiX3 + H/6 x? 3 - 17/6X1X13 - 43/36x2] + n 2 6 [-15/64 XvXr} 

- 59/384 XlX3 + 65/384 x? 3 - 1/2X1X13 - 43/128x2] - 48L^ X iXi3 - 72L?x? - 8Lg 2 x? 3 
+ A(x P ) tt 16 [-1/8 xii?^ + 1/16 xi^p - 1/12 R p qnri Xi3 + 1/48 R c p Xq + 1/48 R c p X i 3 ] 

+ A(x P ) ^ [8/3 ^ x P + 2/3 R c p Xp + 2/3 i$ + A( Xp ) L\ [2/3 i^, Xp + 5/3 R% Xp + 5/3 R d p ] 

- A( Xp )Ll [2 Xl RP qjrv - Xl R c p ] - A( Xp )L r 5 [2/3 R% v Xl3 + 1/3 ^ X is] + ^(x P ) 2 [1/16 + 1/72 (i?^) 2 

- 1/72 + 1/288 (i?p 2 ] - A( Xp )A( Xps ) [1/36^ + 5/72^-7/144^] 
-i( X p)A(x 95 ) [1/36^ + 1/24^ + 1/48^] - i(x P M(Xm) [l/72 R p qmn R v mnl3 - 1/144 ^„ 13 ] 
+ l/8A( X p)A( Xl3 ) - 1/108 A{ Xp )A{ Xst )(Wi mi ) 2 R z tsspp + A( Xp )B( Xp , Xp ,0) [1/4 Xp - 1/18 R p qnT ,R c pXp 

- 1/72 i?^^ + 1/18 (i?p 2 xp + 1/144 R c p R d ] + 1/18 A( Xp )B( Xp , Xm , 0) R z qmp R c p R v mnl3 Xp 

- Mx P )B( Xq , Xq ,0) [l/72R p qnr] R d q -l/U4R c p R d ] - 1/12 ^( Xp )B(Xp S , Xp S , 0) R^ X P s 

- l/lSA( Xp )B(xi, X3 ,0)R g pwri R c p Xp + l/lSA( Xp )C( Xp , Xp , Xp ,0)R c p R d pXp + A{ Xp ,e)ir w [1/8 

- 1/16 Xl R c p - 1/16 R c pXp - 1/16 i^] + A( Xps )tt 16 [1/16 Xps - 3/16 Xgs - 3/16 X i] - 2A( Xps )L£ Xps 

- 5A(xp S )LSXpa - 3A( Xps )i^ Xl + A( Xps )^Xi3 - l/32A( Xps ) 2 - 1/16 A( Xps )A( Xpt ) 

+ l/16A( Xps )A( X9t ) + A( Xps )A( Xm ) [7/144 R™ pp - 5/72 R™ s - 1/48 R™ q + 5/72 i?™s - l/36iZ£ 3 ] 
+ i(Xp S )B(Xp,X P ,0) [l/24fl? w ,xp- 5/24 J^XpJ - A( Xps )B( Xp , Xm , 0) [l/18 #™ s ^ pm Xp 
+ 1/9 R™ s R z qpmXps ] - l/48A( Xps )B( X9 , X9 ,0)i?^ + 1/18 A{ Xps )B{ Xu X3 , 0) R% Xs 
+ l/9A( Xps )B( XuX3 ,0,k)Ri smi + 3/16A( Xps ,e)ir 16 [ Xs + Xl ] + A( Xm )^ 16 [l/16 Xl iC„i 3 

- 1/48 iCm 3 Xm + 1/16 iCmsXis] + A( Xm )LS [4i?™ 13Xm + 2/3i?^„ 1 3 Xm ] - 8A( Xm )L^ Xm 

- 2A( Xm )^ Xm + A( Xm )L r 3 [4i?™ 3Xm + 5/3i?^„ 13Xm ] + A( Xm )i^ [4 Xm + xi^m„i 3 ] 

- A{ Xm )Ll [4/3K"i 3 Xi3 + l/3^„ 13 Xi 3 ] + 1/288 A( Xm ) 2 (iC nl3 ) 2 - 1/108 A( Xm )A( Xst ) R v mnl3 R v mns 
+ A( Xm )B( Xp , Xp ,0) [-1/18 RgmnR^pp Xp - 1/108 (R p mn f R^s Xs + 1/18 R™ p R c p Xp + 1/144 R d p R v mnl3 ] 
+ A(Xm)B(Xp,Xm,0) [-l/54-RJ mn i?^, s i?™ s Xs + l/5AR p smn R™ qs R™ ss Xs + 1/18 (i?£p p ) 2 Xp 

- 1/18 -R™pp^™l 3 Xp] + ^(Xm)-S(Xp: Xn, 0) [-1/54 R P mn Rnss^-mps Xs + 1/54 R P m n^nss^-mqs Xs 

+ 1/18 R™ pp R n mpp X P - 1/18 R n n ppR7ni 3 Xp] - 1/12 A( Xm )B( Xps , Xps , 0) fl™ s Xps 

- 1/108 A( Xm )B( Xm , Xm ,0)(RZ s ) 2 R v mnl3 Xs - l/MA( Xm )B( Xm , X n,0) [R™ s KL^C ps X, 

- RnqsRnss^mps Xs] — 1/108 -A(Xm)-B(Xn, Xm 0) ^nss^mss-^nml3 Xs 

- -4(Xm)-B(xi, X3, 0) [1/18 i?^ m „i?™ p Xp - l/54-R™ s i?s mn i?f m „ Xs] 

+ l/18i(x m )C(xp,Xp,Xp,0)i 1 , ™ p ^Xp + A( Xm;£ )7r 16 [l/8x ro -l/16xi-Rmni 3 -l/8^i 3 Xm 

- l/16iC„i 3 X m ] + l/l^A{ Xv )A{ Xv )Rl nls R^ l3 + A( Xl )A( X3 ) [-1/72^^ + 1/36^^ 
+ 1/144^] - 4A( Xl3 )iI X i 3 - 10A(xi 3 )^Xi3 + l/8A( Xl3 ) 2 - 1/2 A( Xl3 )B( Xl , X3 , 0, k) 

+ l/4A(xi 3 ,e)_7ri 6 xi 3 + 1/16 A( Xls )l( X3s ) - 16A( Xst )L^ Xst - 4A( Xst )i2Xst + 8 A( Xst ) L r 4 Xst 
+ l/18A( Xst )B( Xp , Xp ,0)RP mj Rl nXst + l/54i( Xst )S( Xp , Xm ,0) [R^ t - R^ t ] 

- l/54A( Xst )B( Xp , Xm ,0,k) [R^t-R^t] + l/18A( Xst )B( Xm , Xm , 

0) RnstRmnl 3 Xst 

+ l/54A( Xst )B( X7r , Xn ,0) [R;P-RZ1] - l/54A(x, t )B( Xw ,x„0,fe) [^ - K^t] 

- l/54A( Xs t)B(xi,X3,0)^r 3si + l/54A( Xst )S( Xl , X3 ,0,fc)i?5 t 3st + l/4% sf , e )^ 16Xsf 
+ S( Xp , X p, 0) tt 16 [1/16 X i^ + 1/96 R d Xp + 1/32 R d Xq ] + 2/3 S( Xp , X P , 0) ^p Xp 

+ 5/3 B{ Xp , Xp , 0) LI R d Xp + B( Xp , X p, 0) [-4 Xl i?^ Xp + 2 Xl i?^ Xp + 3 Xi-Rp] 

+ B(X P , Xp, 0) L r 5 [-4/3 J^, xiX3 + 2/3 R c p X \ - 1/3 R d Xl3 ] + B( Xp , Xp , 0) L r 6 [8 X ii?^„ Xp - 4 Xl i? p Xp 

- 4 Xl i$ +4 B( Xp , Xp , 0) U 7 (R d ) 2 + B( Xp , X p, 0) [8/3 X iX 3 - 4/3 R c p X 2 p + 4/3 (R d ) 2 ] 
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+ B(X P , X P , O) 2 [-1/18 R p q ^R d p X P + 1/18 R c p R d x P + 1/288 (R d p ) 2 ] 

+ l/l8B( Xp ,Xp,0)B(Xp,Xm,0)R z qmp R d p R v mni:i Xp ~ 1/18 B( Xp ,x P ,0)B(xuX3,0) R q pnr ,R d pXp 
+ l/l8B(x P ,Xp,0)C(xp,Xp,Xp,0)(R d p ) 2 X P ~ l/16B(x„,X P ,0,e)7ri 6 [xiR d + R d X P ] 

+ 8 B(X P , Xm, 0) Rq pm R p Rm456np + 8 /3 B(Xp, Xm, 0) L\ Rq pm R p Rm456np ~ 6 B(x P s, Xps, °) L\ %1 Xps 

- 2B(x P s,Xps,0)L r 5 x 2 ps + 12 B( Xps , Xps, 0)L r eX i Xps + 4 B( Xps , X P s, 0) L r 8 X 2 ps 

+ 2B(X mi Xmi Q)LlxiR V mnlzXm + 2/3 B( X mi Xmi m, Xm, 0) Lq, XlRmnl3 Xm 

+ 4B(Xm,Xm,0)£ 7 ^ 31 i mm (-R m 45 6n i) + 4/3 B(Xm , Xm , 0) L\ [-R 3 H TOTO (-R m 456nl) ~ ^™13 Xm] 
+ 8 B(x-K,Xm 0) L T 7 ^3117^^456^1-^45671-1 + 8/3 £?(Xtt, Xr,, 0) £g ^3117^-^4567/1 -^456tt1 

+ 1/144 B( Xl , Xl ,0)B(x3,X3,0) J Rfi 1 , 3 d - SB(xi,X3,0)L r 7 RiR d - 8/3S( X i,X3,0)Lg^ 

7 Xp ! Xp J 0)^Xi-RpXp + 4/3C( X p , Xp 7 Xp 7 0)L r 5 R d pX 2 p - 8C( Xp 7 Xp 7 Xp 7 

0) L\x\R d p Xp 

- 8/'SC(x P ,Xp,X P ,0)L r s R d x 2 p + ff F (l,Xp,Xp,Xi3,Xi 3 ) [1/8 Xp - 1/16 Xis - 1/36 (i^) 2 xi 3 

+ 1/36 i?^i? p xi3 - 1/144 (i?p 2 xi 3 ] + H F (1, x P , Xis,X3s,Xu) [-1/12 R p q7rv X q s + 1/24 RP qnr) Xw 

- 1/16 R p swr) Xp + l/A8R p svr)Xq + 1/24 R c pXp s] - l/8H F {l,x P s,X q uXsuX^)Xst 
+ H F '(l,Xm,Xp7Xi3,Xi3) [1/36 -Rg m „-Rm„i 3 Xi3 - 1/72 i?p-Rm„i 3 X13] 

- 1/144 ^(l,Xm,X m ,Xl3,Xl3)(-R™„i 3 ) 2 Xl3 + 1/48 /ffl^m, Xis, X3s,Xl 3 ) [-RnpsXp + RnpsX q 
+ R-mnps Xm + R m nl 3 Xs] - 1/72 H F (1, Xtt, Xr;, Xl3, X13) ^7/13-^71-13 Xl3 

+ H F {1, xi, Xi3, X3, X13) [1/36 RP qn7] R c q X13 - 1/18 R\^R\^ Xi3 - 1/72 R\Rl X is] 

- l/8ff F (l,Xi3,Xi3,Xi3,Xi3)Xi3 + ff F (2,X P ,Xp,Xi3,Xi3) [1/36^^X13-1/72^^X13] 

- l/72/f F (2,x P ,Xm,Xi3 , Xl3) -Rp-Rm„l3 Xl3 + # F (2, Xp,Xl3,X«,Xl3) [1/36 R^rjRp Xl3 - 1/72 i?gi?pXl3] 

+ 1/24 ii^(2,Xp,Xis,X3 S ,Xi3)i?p Xps - 1/144 ff f (5,Xp,Xp,Xi3,Xi3)(-Rp) 2 Xi3 

- l/72 J ff F (5, X i,X3,Xi3,Xi3)^i?3Xi3 + ^ F '(1,Xp,Xp,Xi3,Xi 3 ) [-l/8x P Xi3-3/16x 2 3 

- 1/36 (i?^) 2 Xi 3 + 1/12 R p qvri R c p x\z- 1/48 (R c p ) 2 x 2 13 \ + H F ' (1, Xp , Xu, X3s, Xis) [-5/24 i?^ X ,Xi3 
+ 1/12 i?^ Xps Xi3 + 1/16 R p ^ Xq Xi3 + 7/48 R^ X i 3 Xs - 1/24 R c p Xps Xi 3 ] 

+ 1/8 H F '(1, Xps, X q t,Xst,Xi3)xisXst + ff F '(l,Xm,Xp,Xi3,Xi3) [l/9i?^ m „i?™i 3 X 2 3 

- 1/36 i?P m „i? 

mnl3 

X 2 3 + 1/18 f?" p R c p xh + 1/72 f?£f? 

Xis] 

+ if F ' (1, Xm , Xm , Xi3, Xi3) [1/9 (iC 13 ) 2 X 2 3 + V9 ^Cis^Crns X 2 3 + 5/144 (R v mnl3 ) 2 xl 3 ] 

+ H F ' (1, X m, Xis, X3s, Xis) [-1/16 i?" s XpXis - 5/48 i?™ s x g Xi3 + 1/6 R£ 13 X \ 3 - 1/48 fC„ P s x™Xi3 

- l/48iC n 13Xl3X-] + ^'(l,X7T,Xr,,Xl3,Xl3) [1/18^^X13 + 1/9^13^13X13 

+ l/72i?^ 13 ^7ri3:X 2 3] + ff F '(l,Xi,Xi3,X3,Xi 3 ) [-1/36 R^R c q x? 3 + 1/6 R\^Rl vri x\ 3 

+ l/72f?^ X 2 3 ] + l/8ff F '(l,Xi3,Xi3,Xi3,Xi3)x 2 3 - # F '(2,Xp,Xp,Xi3,Xi3) [l/36i?^Xx 2 3 

- 5/72i?^x? 3 ] + # F '(2,x P ,X m ,Xi3,Xi3) [1/18 K Pp R d pX% + 1/72 R d R v mnl3 xl 3 ] 

- H F '(2, Xp ,Xi3,X q ,Xi3) [1/36 Rl mi R d p X 2 13 - 1/72 R c q R d x 2 13 ] - l/24H F ' (2, Xp ,Xis,X3s,Xi3) R d p X P sXi3 
+ 5/U4H F ' (5, X p,Xp,Xi3,Xi3) (R d ) 2 X 2 i 3 + l/72H F '(5, X i,X3,Xi3,Xi3)R d R d 3X 2 i3 

+ iff (1, x P , X P , Xis, Xis) [x 2 3 + 2/9 (R p qwv ) 2 xi 3 - 2/9 R p ^R c p x? 3 + 1/9 (R c p ) 2 xi 3 } 

+ 1/3 H[ (I, Xp, Xis, X3s,Xl3) R q -x n R Z s q p Xl3 - 1/3-fff (l,Xps,X<7SiXp,Xl3)-RLr,Xl3 

- 1/3 iff (1, Xps, X«s,Xm,Xl3) ^13-^^X13 + 1/9-fff (l,Xl3,Xp,Xm,Xl3) [-R^ mn -Rmnl3 Xl3 

- -^rTpp -Rgpm -Rp X13] - (1) Xl3, Xm, Xm, Xu) [l/9 -Rr?13-^mnl3 Xl3 + 1/18 (-Rmnl3) 2 Xl3] 

- 1/9 iff (1, Xl3, Xtt, Xr;, Xl3) [R^ppR^pR^TrlS Xl3 ~ ^13^x13 X13] 

+ Fr(l,Xi3,Xi,X3,Xi3) [l/9RP q ^R c q x 2 13 -2/9Rl 11 Rf 7rri x 2 i 3 ] 
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+ 1/9 iff (3,xi3,X P ,Xp,Xi3) K m) ^X?3-^X?3] + l/9H?'(3,xi3,X P ,X q ,Xi3)R% nr ,RU 2 i3 
+ 1/9 iff (3, Xl3,Xp,Xm,Xl3)iCl 3 ^ gm ^pXl3 - l/18fff'(7,Xi3,Xp,Xp,Xi3)(^) 2 X?3 

- 3/8 i?^ (!,Xp)Xp)Xl3,Xl3)Xl3 ~ 1/ 8 #21 ( 1 :Xp I Xls,X3s,Xl3)-R9^-Rf 9p Xl3 

+ l/8H2i(l,XpB>Xqa,Xp>Xl3) [Rqirq Xl3 + R^nr, Xl3 _ -^p X13] + V 8 #21 (!) Xps, Xgs, Xg, Xl3) R^Rlpq Xl3 
+ 1/ 8 #21 (l'XpS)X«s,Xm,Xl3)-R™13-Rpgm^spmXl3 _ 1/ 8 #21 (l) Xm, Xls> X3s, X13) -R™13^slm-^s3m Xl3 

+ #Ti (1, Xi3, X P , Xp, X13) [3/16 X 2 3 + 1/12 (^) 2 Xi 3 - 1/12 R p qmi R c p X 2 3 + 1/48 (f? P ) 2 X 2 3 ] 

- #21 (!>Xl3,Xp,Xm,Xl3) [l/12^ mn i?m„i3Xl3 ~ l/24-Rpi?m„i3Xl3] 

+ 1/48^(1 

, Xl3, Xm, Xm, Xl3) (^mnl3) 2 Xl3 + 1/24 #21 (1, Xl3, Xtt, Xr,, X13) ^,,13-^^13 Xl3 

+ Hi! (1, X13, Xi, X3, X13) [-1/12 R p q7IV R c q xl 3 + 1/6 Rl^RU X 2 3 + 1/24 R C X R C 3 x\ 3 ] 
+ 3/8iJ 2 p 1 '(l,xi 3 ,xi3,xi3,xi 3 )x 2 3 + 3/8^ (l,x s t,Xp S ,X<jt,Xi3)Xi3 

- l/8ff£'(3, Xps , X p,X 9S ,Xi3)i?pX 2 3 - #21' (3, Xi3 ,Xp,Xp,Xi 3 ) [l/12 J^Jj£ x? 3 - 1/24 R c p R d p X ? 3 ] 

- H^'(3 ) Xl3,Xp.Xg,Xl3) [l/12^ ff Xx 2 3-V24^Xi 3 ] 

+ l/24i/ 2 i? 1 '(3,xi3,Xp,Xm,Xi3)i?X-i3Xi3 + l/48^(7,xi3,Xi»X P ,Xi3)(i^) a Xi3 

+ l/24i/ 2 F 1 '(7,xi3,Xi,X3,Xi3)^i?3 d X 2 3- (73) 




FIG. 3: The relative shifts of the charged pseudo-scalar meson mass Am and decay constant Af to NNLO for d va i = 1 and 
dsca = 1, as a function of the valence and sea-quark masses xi an d X4- The quantity plotted represents the sum of the NLO 
and NNLO shifts, and the difference between two successive contour lines in the plots is 0.10. The values chosen for the LEC:s 
correspond to "fit 10" as discussed in the text. 



V. NUMERICAL RESULTS AND DISCUSSION 

A. Checks on the Calculation 

It is mandatory, in such a lengthy and complex cal- 
culation, that all possible measures be taken to ensure 
the correctness of the analytical as well as the numeri- 



cal results. First of all, the divergence structure and the 
cancellation of nonlocal divergencies has been found to 
behave as expected. However, finiteness alone is not a 
very rigorous check on the end results, as several of the 
two-loop sunset integrals are convergent because of the 
appearance of double poles in the propagators. A much 
more rigorous check is provided by the fact that most 
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FIG. 4: Comparison of the NLO, NNLO and total NLO+NNLO shifts of the charged meson mass for d va i = 1 and d aca _ = 1. 
The left-hand plot shows the NLO and NLO+NNLO results for the set of LEC:s labeled 'fit 10', which has been used for 
all the other figures in this paper. The right-hand plot shows the result when "fit 10" is augmented by Lq — —0.2 x 10 -2 , 
LI = 0.1 x 10" 3 , and L r 6 = -0.1 X 10" 3 . The effects of introducing a nonzero value for the 0(p 6 ) LEC TfJs is a l s ° demonstrated 
for the NNLO and NLO+NNLO results. 



parts of the calculation have been performed indepen- 
dently by each one of the authors, such that the computer 
programs written to handle the symbolic manipulations 
(mainly using FORM) have been developed independently. 
At the end of the calculation, and at several intermediate 
steps during the process, the output from each program 
has been compared and cross-checked, to make sure that 
perfect agreement was found. These checks are highly 
nontrivial, since the form of the (equivalent) outputs is 
by no means unique. Firstly, the exact appearance of the 
result depends on the order in which the various sym- 
metries inherent in the Feynman diagrams have been 
implemented. Secondly, the various loop integrals and 
propagator residues R satisfy a large number of nontriv- 
ial relations, yielding further possibilities to rewrite the 
result. In the very end, the output from one computer 
program was simplified and compressed (by more than an 
order of magnitude) to bring it into a publishable size. 
Again, that expression was thoroughly checked to de- 
termine whether the compactified analytical expressions 
were still equivalent to the other two versions. 

The computer programs for the numerical treatment 
of the analytical expressions have also been produced in- 



dependently of each other. Thus the risk for acciden- 
tal cross-contamination has been minimized as each code 
was produced in a different programming language (f 77, 
f90 and C++). All three programs have been found to 
agree with each other, and the permutation symmetries 
under exchange of quark masses have been verified up 
to the numerical precision. As a final check, the general 
mass cases have been found to converge numerically, as 
expected, towards the less complex ones when pairs of 
quark masses become degenerate. 



B. Presentation of Numerical Results 

In view of the complexity of the analytical NNLO 
results in PQxPT, a graphical presentation of the re- 
sults is clearly called for. In general, the corrections 
to the pseudoscalar meson masses and decay constants 
are real- valued functions of all the sea and valence quark 
masses. Since, for the nondegenerate cases, this involves 
the graphical representation of a function which depends 
on up to five different quark masses, an exhaustive plot 
is obviously difficult to produce. Except for the most de- 
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FIG. 5: The combined NLO and NNLO shifts A M of the 
charged pseudoscalar meson mass, plotted for d va i = 2 and 
dsca = 1, for 8 — 60° in the xi ~ Xi plane, and for the pro- 
portionality factors x — {0.25, 1.1, 2.0, 5.0, 10.0} between the 
valence quark masses X3 an d Xi- 



generate (1+1) mass case, the choice has therefore been 
to present everything as a function of the valence quark 
mass Xi on lyi by parameterizing the dependence on all 
other quark masses in terms of xi- Furthermore, due 
to the size and complexity of the analytical expressions, 
their practical usefulness depends highly on the availabil- 
ity of software which can produce numerical output from 
these expressions in a reliable manner. In view of this, 
the software which was used for the figures in this pa- 
per will be made available in the near future from the 
website Q- 

In addition to the quark mass dependence, the NNLO 
expressions are also functions of a number of largely un- 
known LEC:s. In the long run, these LEC:s should of 
course be determined by a fit of the PQxPT formulas to 
Lattice QCD data. At the present time, this is not yet 
possible, although suitable simulation results should be- 
come available in the near future. Therefore, the present 
work makes use of the LEC:s determined by a fit to exper- 
imental data, referred to as "fit 10", which is presented 
in Ref. H3|. That fit has F = 87.7 MeV and a renor- 
malization scale of fi = 770 MeV. The NNLO LEC:s K\ 
and the NLO LEC:s L\, L r e and Lq were not determined 



in "fit 10", and they have thus been set to zero for sim- 
plicity. However, some results for nonzero values of these 
LEC:s are presented in Fig. 0] It should be noted that 
Lq cannot, as discussed earlier in this paper, be deter- 
mined from experimental data, since it is a distinguish- 
able quantity only in the PQ theory. Some recent results 
on L\ and Lg have been obtained in Ref. |2~j|. but they 
have nevertheless been set to zero in most of the plots 
in this paper, since the present numerics are mainly in- 
tended for illustrative purposes. 

In the next subsections, the NNLO meson masses and 
decay constants are presented in terms of the relative 
shifts Am and respectively, which represent the 
change in the indicated quantity due to the NLO and 
NNLO contributions. For the pseudoscalar meson mass, 
Am is defined by 

A M = M p 2 hys /xy - 1, (74) 

where Xij has again been substituted for the lowest order 
result. The calculation of Am thus involves the numerical 
evaluation of Mph ys in Eq. (|52H to the order of the plot 
in question. For clarity, it is useful to recall here the 
definitions 

Xi = 2B m qt , 

Xij = (Xi+Xi)/2, (75) 

where m q i denotes the (current) mass of the quark c^. 
In xPT, this corresponds to the lowest order mass of 
a meson composed of a quark and an antiquark, each 
of mass m q i. As an example, the quantity Xl plotted 
on the horizontal axes of Figs. |3| and 0] represents the 
lowest order mass of the valence quark meson qiqi, and 
is given by xi = 2 Bo 771 1- Where no confusion can arise, 
the quantities Xi are sometimes referred to as "quark 
masses" . 

The decay constant shift Ap is defined in a similar 
way, such that 

Ap = F phys /F - 1, (76) 

where F p h ys of Eq. I|70|) is again evaluated to the desired 
order. Since the decay constant has been treated in de- 
tail in Ref. [Tol ] , the numerical analysis of this paper will 
mainly focus on the pseudoscalar meson masses. The de- 
pendence of the results on the sea-quark mass Xi and 
the valence quark mass X3 is parameterized in terms of 
an angle 6 in the xi ~ Xi plane and a proportionality 
factor x, according to 

Xi = tan 9 xi, 

X3 = xxt, (77) 

and the deviations of the remaining sea quark masses 
from Xi are similarly given by 

xs = yxi, 

Xe = zxi- (78) 
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FIG. 6: The combined NLO and NNLO shifts Am of the charged meson mass plotted for d sca = 2 and 8 = 60°. The left-hand 
plot shows the results for d va i = 1 and a ratio between the sea-quark masses xs an d X4 of z — {0.25, 0.5, 1.0, 1.1, 3.0, 10.0}. 
The right-hand plot shows the d va i = 2 result for z — 3.0 and a ratio between the valence quark masses X3 an d Xi °f 
x = {1.0, 1.1, 2.0, 3.0, 3.5, 5.0}. 



The values of x,y,z and 9 shown in the various plots 
in this paper have been chosen for convenience, and do 
not carry any particular significance. A value 9 > 45° 
was chosen since the sea-quark masses are then heavier 
than the valence ones, a situation which is encountered 
in realistic Lattice QCD simulations. The values for the 
other quark mass ratios have been picked so as to illus- 
trate a large range of possible quark mass combinations, 
but they have otherwise been arbitrarily chosen. 

The expressions calculated to NNLO for degenerate 
quark masses should be numerically recoverable as lim- 
its of the more general ones, which serves as a useful 
consistency check. This has been checked for all cases 
studied, and is explicitly demonstrated for each of the 
proportionality factors x, y and z. For each diagram, one 
of the plot curves from a more degenerate mass case (the 
thicker black curve) has been included for comparison. 
The curve with the new mass ratio, x, y or z equal to 1 
shows the more degenerate case and the curve with the 
ratio equal to 1.1 shows how the more degenerate case is 
approached. 

A naive validity criterion for PQ;\;PT is that all masses 
should satisfy approximately Xi < 0-3 GeV 2 , but this 
should naturally be separately checked for each quantity 



in question. Note that the plots presented in this paper 
have been extended to include also regions where this 
constraint is not satisfied. A general study of the conver- 
gence of the PQxPT expansion up to NNLO is beyond 
the scope of this paper, as realistic values obtained from 
Lattice QCD simulations should be used for all LEC:s 
before any meaningful statements concerning the relative 
magnitude of the NNLO corrections can be made. The 
question of convergence then hinges on whether it is pos- 
sible to describe the Lattice QCD data in such a way that 
the NNLO correction is of a reasonable magnitude. How- 
ever, on a qualitative level, the convergence is typically 
better for the decay constant than for the mass, and the 
inclusion of nonzero values for the L\ tends to improve 
the convergence. Nevertheless, it should be kept in mind 
that these statements are mainly valid for the LEC:s of 
"fit 10" since a general study of the convergence has not 
yet been performed. 

It should also be noted that an example which demon- 
strates the effect of variation of the 0(p 4 ) LEC:s as well 
as the introduction of nonzero values for the 0(p 6 ) LEC:s 
on the convergence of the chiral expansion (in the most 
degenerate mass case) is shown in Fig. 0] and discussed 
in Sect.lVCl 
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FIG. 7: The combined NLO and NNLO shifts of the charged meson mass, in the left-hand plot for d va l = 1 and d aGa = 3, 
with 6 — 60°, z — 3.0 and a ratio Xs/x^ °f V = {0-5, 1.0, 1.1, 2.0, 4.0, 8.0}. The right-hand plot shows the case of d va i = 2 and 
ioa = 3, with 6 = 60°, z = 3.0, x = 2.0 and y = {0.5, 1.0, 1.1, 4.0, 10.0} 



C. Numerical Results for d sca = 1 

Numerical results for the simplest possible case, where 
both valence quarks are degenerate, i.e. x — 1, have 
already appeared in a number of earlier works. In partic- 
ular, the case of degenerate sea-quark masses, y = z = 1, 
has been considered extensively in Refs. and [ljj both 
for the meson mass and the decay constant, as well as 
for different choices of LEC:s. In addition to this, Fig- El 
shows a two-dimensional plot of the mass shift Am and 
the decay constant shift Ap over the whole parameter 
space up to 0.5 GeV 2 . These include thus Am and Ap 
for all values of 9. 

An important question in the context of PQ^PT calcu- 
lations is the issue of convergence of the chiral expansion 
as a function of the input quark masses. At present, 
this question cannot be easily answered since the behav- 
ior of the NNLO expressions given in this paper depends 
rather sensitively on the values of the LEC:s. This point 
is well illustrated by Fig. 21 where the NLO mass shift for 
9 = 60°, d va \ = 1 and rf soa = 1 has been given explicitly. 
For the set of LEC:s labeled "fit 10", the NLO result is 
very small and the total result is completely dominated 
by the NNLO contribution. For this set of LEC:s, there 



is obviously no good convergence and the overall correc- 
tions are rather large for most values of xi- Furthermore, 
the strong curvature exhibited by the NNLO result is not 
typical of the behavior seen in Lattice QCD simulations. 
It should be noted that unquenched x?T also shows sim- 
ilar behavior (23 for these values of the LEC:s. 

This apparently pathological behavior can be com- 
pletely changed by appropriately tuning the values of the 
LEC:s, as shown in Fig.QJ From the results shown in that 
plot, one can conclude that "fit 10" with small changes 
is actually completely compatible with excellent conver- 
gence of the PQxPT expansion. In fact, by choosing 
L 7 = -0.2 x 10~ 2 , L\ = 0.1 X 1Q- 3 and L r 6 = -0.1 x lO" 3 , 
the dramatic curvature in the NNLO contribution is 
straightened out and the near total cancellation in the 
NLO result is eliminated. The total NLO+NNLO re- 
sult then lies close to the NLO one even for rather large 
values of xi arid Xi- The effects of considering nonzero 
values of the LEC:s at 0(p 6 ) can be sizeable as well. The 
changes that result from the introduction of a naturally 
sized value of K^ 5 = —0.3 x 10~ 5 are also shown in Fig.01 
Overall, these changes in the values of the LEC:s yield a 
nicely convergent chiral expansion and a rather smooth 
and featureless final result. 
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FIG. 8: The combined relative NLO and NNLO shifts of 
the charged meson decay constant, Af, plotted for d va i = 2 
and d soa = 3, with 6 = 60°, z = 5.0, x = 2.0 and y = 
{0.25, 0.5, 1.1, 1, 2, 4}. Note that further plots for smaller val- 
ues of d va i and d Bca can be found in Ref. [Tl|. 



X3 is given by the ratio x as denned in Eq. (f77|l . For the 
case of c£ soa = 2, the behavior of Am as a function of xi 
is illustrated, for different choices of x and z, in Fig. 
As for the previous mass case, all mass shifts vanish in 
the chiral limit for all cases shown in Fig. H3 In the left- 
hand plot of Fig. the thick black curve of the mass 
case (1+1) has been plotted for the same parameters as 
in Fig. and as expected, the NNLO correction again 
approaches this curve in the limit z — > 1. In the right- 
hand plot of Fig. El the thick black curve corresponds to 
the mass case (1+2), with 60° and z = 3.0. There, the 
convergence to the more degenerate quark mass configu- 
ration is quite fast, as the plot for x = 1.1 lies very close 
to that curve. 

Finally, the case with d sea = 3 depends on all three 
sea-quark mass parameters X4j X5 an( i X(n given in terms 
of the 6 and y,z defined in Eqs. l(77|) and lf7%|) . In the 
left-hand plot of Fig. [7| Am is shown as a function of 
Xi for (i V ai = 1 and d sca — 3, and in the right-hand plot 
the case with d va \ = 2 and <i S ca = 3 is shown. As for 
all previous cases, the mass shifts vanish in the chiral 
limit. In both plots of Fig. [7| the thick black curves 
can be found among the lines plotted for d soa = 2 in 
Fig. |SJ Again, the convergence to the respective thick 
black lines is seen to be properly attained when y — ► 1. 
For completeness, a plot of the relative NNLO correction 
to the decay constant in shown in Fig. [5] Plots for the 
more degenerate mass cases can be found in Ref. |10). 



E. Determination of LECs from Lattice QCD 
simulations 



For the case of d va i = 2 and d sca = 1, the mass shift 
Am has been plotted in Fig. [S] as a function of the input 
valence quark mass xi fo r several values of the valence 
quark mass ratio x. Note that the shifts A should vanish 
in the chiral limit, Xi ™> 0, which is obviously satisfied 
by the plots since the sea-quark masses also vanish for 
Xi — > 0. On the other hand, it is apparent from Fig. |3| 
that the shifts A do not vanish if the valence quark mass 
is set to zero for a constant value of the sea-quark mass. 
This effect is due to the quenched chiral logarithms and 
constitutes a generic feature of the PQ theory. Further- 
more, Fig. [S] shows that the results for the (2+1) mass 
case converge well toward the (1+1) mass case (the thick 
black curve) in the limit x — > 1. 



D. Numerical Results for <i sca = 2 and d sca = 3 

The expressions for d se a = 2 have two sea-quark mass 
parameters, \i an d Xe- l n the numerical analysis, the 
latter is given in terms of the former by the ratio z of 
Eq. (|75|) . As for d sca — 1, the sea-quark masses are pa- 
rameterized with respect to the valence quark mass xi 
by the angle 9, as in Eq. (|77fl . Similarly, for d va \ = 2 the 
relationship between the valence quark masses \i an d 



As noted in the previous sections, the main motivation 
for NNLO calculations in PQxPT is that such expres- 
sions can be used to analyze the results of PQ Lattice 
QCD simulations in terms of the LEC:s of unqucnchcd 
QCD. Up to this point, the numerical results presented 
in this paper have illustrated the quark mass dependence 
of the NNLO expressions for a given set of LEC:s. In 
this section, the attention is turned towards finding effi- 
cient methods of determining the LEC:s from fits of the 
NNLO expressions to future Lattice QCD data. It should 
be kept in mind that the LEC:s denoted by L\ and K\ 

always refer to the L^ 3p9 ' and the K^ 3pqS> of the PQ the- 
ory with n soa = 3. An equivalent discussion for the case 
of risen = 2 can be found in Ref. 11]. 

At NLO, the expressions for the masses and decay con- 
stants of the pseudoscalar mesons can be separated into a 
tree-level contribution, which is a function of the L^, and 
a remaining part independent of the L\, which involves 
the chiral logarithms. In the following developments, the 
average quark mass is denoted by x, where the index de- 
noting the power of the sea quark masses averaged has 
been dropped for simplicity. Note also the use of the 
product of lowest order sea-meson masses x^Xn^ wmcn 
was defined in Eq. 1)28)1 . For the most general case of 
e? va i = 2 and d sca — 3, the dependence on the is then 
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proportional to 

/ c ( t 4)23 ~ 3xL r 4 + xi 3 L r 5 (79) 

for the decay constant, and 

<^ )23 ~ 3x (2L r e - L\) + X i3 (2L r s - LI) (80) 

for the mass 0, from which an overall factor of Xi3 nas 
been removed. The expressions for the more degenerate 
quark mass cases can then be obtained by consideration 
of the appropriate limits. For the expressions in Eqs. I|79|l 
and (|80() , these limits can be taken straightforwardly, and 
it can be seen at once that only the simplest (1+1) mass 
case is required in order to determine all the NLO LEC:s 
present in the expressions, since one can determine L\ 
and L\ by fitting Eq. (|79|) to Lattice QCD data, and 
then use this knowledge to determine the remaining two 
constants from Eq. I]80p. The more complicated cases 
with nondegenerate quarks are thus redundant for a de- 
termination of the LEC:s at NLO. 

A similar analysis at NNLO becomes more challenging, 
not because of the analytical complexity of the NNLO ex- 
pressions, but rather because the chiral logarithms are no 
longer independent of the L\. On the other hand, the de- 
pendence of the NNLO expressions on the K\ is similar 
to that of the NLO ones on the L\, with the exception 
that the K\ parameters at NNLO are slightly more nu- 
merous than the L\ at NLO. In all, the 0(p 6 ) expressions 
depend on a total of 12 K\ parameters, all of which occur 
in the expression for the meson mass, whereas the expres- 
sion for the decay constants only depends on five of them. 
For the decay constants, the expression for d va i = 2 and 
^sea — 3 can be factorized into the form 

/ c ( t 6)23 ~ 2 X i 3 K r 19 - X iX3(K r 19 -K r 23 ) 

+ 3xXi3K r 20 + 9 X 2 (K r 21 +K r 22 ) 
-Qx,X n K r 21 , (81) 

up to an overall numerical factor which has been omitted. 
From this expression, it is then possible to determine all 
5 LEC:s K\ 9 through K 23 . The (1+1) mass case suffices 
to determine three combinations of LEC:s from Eq. I]8ip. 
but in order to separate K[ 9 from K 23 at least one case 
with d va i = 2 has to be considered, otherwise only the 
combination K[ 9 + K 23 is accessible, as can be readily 
seen by setting X3 ~~ * Xi m Eq. I|81l) . Similarly, in order 
to separate K 21 from K 22 , at least one case with d sca = 2 
has to be considered, since only the combination K 2l + 
3K 22 can be distinguished for completely degenerate sea- 
quarks, i.e. d sc& = 1. This analysis is identical in form to 
the one for PQxPT with n sca = 2 given in Ref. [Tj|. It 
is actually possible, by making the replacements Xl3 ~ * 
X12, X -> 2/3x34 and x^Xn -> X3X4/3 to recover the 
corresponding expression for n sca = 2 from the present 
one given in Eq. I|81l) . 

The dependence of the meson masses on the K\ is 
somewhat more involved, but the structure in the quark 



mass combinations is similar to that of Eq. I|81|l . Ex- 
plicitly, the dependence of the NNLO shift of the pseu- 
doscalar meson mass on the K\ is proportional to 



<^ )23 ~ - 2 xl 3 (Klj + K{ 9 - 3K r 25 - K r 39 ) 
+ XiX3 (K[ 9 - K r 23 - 3K r 25 ) 

- 6 x Xi3 (K r ls + K r 20 /2 - K r 26 - K$ ) 

- 9 X 2 (K r 21 +K r 22 -K r 26 -3K r 27 ) 

+ Qx.Xv (K r 21 ~K r 26 ), (82) 

from which a factor of Xi3 as well as an overall numerical 
factor has been removed. As for Eq. (JSTJ, the (1+1) 
mass case here gives access to 3 combinations of LEC:s, 
plus an additional combination for each of the (1+2) and 
(2+1) mass cases. In all, 5 additional combinations may 
be determined from the mass expression, even though the 
dependencies are slightly more entangled. However, it is 
clear that simulations with nondegenerate quark masses 
are needed in order to obtain independent information 
about as many of the K\ as possible. On the other hand, 
when one determines the L\ from NNLO fits, it may be 
convenient to work in the (1+1) mass case only, since it 
is then still possible to distinguish all of the L\, while the 
number of distinct combinations of K\ parameters is as 
low as possible, which makes the fitting somewhat less 
complicated. 

It is also useful to compare the three-flavor result in 
Eq. I|82|) with the analogous formula in the two-flavor 
treatment of Ref. 11]. As for the 0(p 6 ) terms from 
the three-flavor decay constant, the terms from the mass 
expressions can be similarly translated into the expres- 
sion for the two-flavor case using the replacements men- 
tioned in connection with Eq. (|81|l . Any apparent dif- 
ferences in the numerical prefactors are then only due to 
the slightly different notational conventions between the 
present work and that of Ref. ^lj. Some additional L\ 
show up in the NNLO expressions, namely L 9 through 
L 3 and It has been noted in Ref. that L 7 can 
be determined from the properties of the double pole of 
the PQ neutral meson propagators. This quantity has 
not yet been calculated up to NNLO. The other four L\ 
are those which are relevant for meson-meson scattering. 
In the present analysis, their values have been used as 
input, but they could in principle be separated from the 
K\ since the terms where they are present depend non- 
analytically on the quark masses. 

While the number of free parameters in the NNLO ex- 
pressions which need to be fitted to the Lattice QCD 
data is certainly large, the above analysis has demon- 
strated that their number is neither unmanageable nor 
overwhelming. In particular, the possibility of producing 
Lattice QCD data for different combinations of sea and 
valence quark masses should simplify matters consider- 
ably. 
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VI. CONCLUSIONS 

In this paper, a complete calculation of the mass of 
a charged pseudoscalar meson to NNLO in PQ^PT has 
been performed, and explicit analytical formulas of a rea- 
sonable length have been given for all relevant degrees 
of degeneracy in the input quark masses. In addition, 
the hitherto missing NNLO result with three nondegen- 
erate sea-quarks for the decay constant of a charged pseu- 
doscalar meson has been provided here. These expres- 
sions could be brought into a manageable size, since the 
residues of the PQ propagators in the neutral meson sec- 
tor satisfy a large number of relations between sums of 
products of ratios of differences of quark masses. As elab- 
orated in Sect. IV Al the results have passed a variety of 
nontrivial consistency checks, analytical as well as nu- 
merical. 

In the numerical analysis of the NNLO expressions, 
various plots for different choices of the input quark 
masses were presented to indicate the typical size of 
the corrections. From these plots, it is certainly evi- 
dent that the total NLO + NNLO corrections may be 
uncomfortably large, unless the input quark masses as- 
sume very small values. However, it is known that the 
NNLO corrections in unquenched xPT are rather large 
as well [2(3, , which is especially true for the case of 
the meson mass. The fact that the partially quenched as 
well as the unquenched expressions depend sensitively on 
the values of the L£ as well as the K\ suggests that this 



problem may be entirely due to the omission of the K\ 
from the numerical analysis. Thus, once a realistic set 
of values for the LEC:s is available, the convergence of 
the expansion up to NNLO should be dramatically im- 
proved. Nevertheless, one may readily conclude that the 
NNLO effects are bound to be nonneglegible for values of 
the quark masses that are presently used in Lattice QCD 
simulations. 

As yet, no fits to available Lattice QCD data have been 
attempted, as such a comparison is beyond the scope of 
this article. It may be noted that an extrapolation of the 
simulation results to zero lattice spacing and to infinite 
volume is necessary for a rigorous comparison with the 
PQxPT expressions, but on the other hand it is possible 
to simply include such effects as extra uncertainties in 
the determined values of the LEC:s. 
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